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d Abstract 

The enormous experimental progress in atomic, molecular and optical (AMO) physics dur- 
ing the last decades allows us nowadays to isolate single, a few or even many-body ensem- 
bles of microscopic particles, and to manipulate their quantum properties at a level of preci- 
sion, which still seemed unthinkable some years ago. This versatile set of tools has enabled 
the development of the well-established concept of engineering of many -body Hamiltoni- 
ans in various physical platforms. These available tools, however, can also be harnessed 
to extend the scenario of Hamiltonian engineering to a more general Liouvillian setting, 
which in addition to coherent dynamics also includes controlled dissipation in many-body 
quantum systems. Here, we review recent theoretical and experimental progress in different 
directions along these lines, with a particular focus on physical realizations with systems 
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of atoms and ions. This comprises digital quantum simulations in a general open system 
setting, as well as engineering and understanding new classes of systems far away from 
thermodynamic equilibrium. In the context of digital quantum simulation, we first outline 
the basic concepts and illustrate them on the basis of a recent experiment with trapped ions. 
We also discuss theoretical work proposing an intrinsically scalable simulation architecture 
for spin models with high-order interactions such as Kitaev's toric code, based on Rydberg 
atoms stored in optical lattices. We then turn to the digital simulation of dissipative many- 
body dynamics, pointing out a route for the general quantum state preparation in complex 
spin models, and discuss a recent experiment demonstrating the basic building blocks of 
a full-fledged open system quantum simulator. In view of creating novel classes of out-of- 
equilibrium systems, we focus on ultracold atoms. We point out how quantum mechan- 
ical long range order can be established via engineered dissipation, and present genuine 
many-body aspects of this setting: In the context of bosons, we discuss dynamical phase 
transitions resulting from competing Hamiltonian and dissipative dynamics. In the context 
of fermions, we present a purely dissipative pairing mechanism, and show how this could 
pave the way for the quantum simulation of the Fermi-Hubbard model. We also propose 
and analyze the key properties of dissipatively targeted topological phases of matter. 

Key words: Open quantum systems, quantum simulation, atomic physics, trapped ions, 
quantum phase transitions, unconventional pairing mechanisms, topological phases of 
matter. 
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1 Introduction 



The extraordinary experimental progress in AMO physics experienced during the 
last decades allows us nowadays to isolate one or few microscopic particles, or 
even many-body ensembles of them, and to manipulate, control and detect their 
quantum states almost perfectly. Harnessing the available tools offers unique possi- 
bilities to extend the customary idea of Hamiltonian engineering to a more general 
scenario, where coherent and controlled driven-dissipative dynamics appear on an 
equal footing. This program comprises different directions. On the one hand, the 
ability to control both coherent and dissipative dynamics constitutes a complete set 
of tools for general open-system quantum simulation, very much in the spirit of a 
truly universal simulator device. On the other hand, the possibility to combine co- 
herent and dissipative dynamics opens the door to novel classes of artificial out-of- 
equilibrium many-body systems without immediate counterpart in condensed mat- 
ter. In this work, we review theoretical and experimental progress on the quantum 
simulation and open-system dynamics of many-particle systems with cold atoms 
and trapped ions from various perspectives. 

Simulation of quantum physics on classical computers is in many cases hindered by 
the intrinsic complexity of many-particle quantum systems, for which the compu- 
tational effort scales exponentially with the number of particles. Thus Feynman's 
vision was to build a controllable quantum device which can be programmed to 
act as a quantum simulator for any quantum system, and would allow one to study 
complex quantum systems, intractable on classical computers, from a wide plethora 
of research fields (Buluta and Nori, 2009). Such a device can be built as an analog 
or digital quantum simulator, and its time evolution can represent a Hamiltonian 
closed system or open system dynamics. In analog quantum simulation one 'builds 
the Hamiltonian directly' by 'always-on' tunable external control fields. Familiar 
examples are cold atoms in optical lattices as analog simulators of Bose and Fermi 
Hubbard models (Lewenstein et al., 2006; Bloch et al., 2008; Dalibard et al., 2011) 
or Rydberg atoms (Saffman et al., 2010) or trapped ions (Schneider et al., 2012; 
Johanning et al., 2009) for the simulation of spin systems. In contrast, in digital 
quantum simulation the initial state of the quantum system is encoded in a regis- 
ter of qubits. For any many-body quantum system with few-particle interactions, 
the time evolution can then be efficiently approximated (Lloyd, 1996) according 
to a Trotter decomposition in small, finite time steps, realized by a stroboscopic 
sequence of quantum gates, as familiar from quantum computing. As we will dis- 



cuss below, the digital simulation approach can be applied to realize coherent as 
well as dissipative many-body dynamics, in particular of open many-particle sys- 
tems involving n-body interactions and constraints, as they naturally appear, e.g., 
in complex condensed matter models, quantum chemistry, high energy physics and 
many-body spin models of interest in the field of topological quantum information. 

In the context of engineering open many-body systems, cold atomic gases offer 
a natural and versatile platform. A large part of current research in this field fo- 
cuses on tailoring specific Hamiltonians, made possible by the precise control of 
microscopic system parameters via external fields. The resulting systems are well 
described as closed quantum systems, isolated from the environment, and rest in 
thermodynamic equilibrium - in close analogy to condensed matter systems. In 
contrast, here we will be interested in a scenario where many-body ensembles are 
properly viewed as open quantum systems, much in the spirit of the setting of quan- 
tum optics and without direct condensed matter analog: A system of interest is cou- 
pled to an environment in a controlled way, and is additionally driven by external 
coherent fields. As anticipated above, via such reservoir engineering driven dissipa- 
tion may then not only occur as a perturbation, but rather as the dominant resource 
of many -body dynamics. In particular, we point out that, while dissipation is usually 
seen as an adversary to subtle quantum mechanical correlations, in proper combi- 
nation with coherent drive, it can act in exactly the opposite way - even creating 
quantum mechanical order. More generally, the results presented below pinpoint 
the fact that the far-from-equilibrium stationary states of such driven-dissipative 
ensembles offer a variety of novel many-body aspects and phenomena. 

Under rather general circumstances, discussed and justified below, the dynamics of 
the many-particle quantum systems we are interested in here can be described by 
the following master equation: ^ 

d,p = -i{H,p]+L{p) (1) 

for the density operator p(0 of the many-body system (Gardiner and ZoUer, 1999). 
The coherent part of the dynamics is described by a Hamiltonian H = Y,aHa, where 
Ha act on a quasi-local subset of particles. Dissipative time evolution is described 
by the Liouvillian part of the master equation, 

£(p) = J]^i(2c,pcl-clc,p-pclc,), (2) 

where the individual terms are of Lindblad form (Wiseman and Milbum, 2009) 
and are determined by quantum jump operators c^ , acting on single particles or on 
subsets of particles, and by the respective rates ja at which these jump processes 
occur. 



Throughout tfiis article we set /i = 1 . 
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While there have been several comprehensive recent reviews on quantum many- 
body physics, quantum simulation and computation with quantum optical systems 
(Jane et al., 2003; Ladd et al., 2010; Cirac and ZoUer, 2012) involving atoms (Dal- 
ibard et al., 201 1; Bloch et al., 2008; Lewenstein et al., 2006; Baranov et al., 2012; 
Bloch et al., 2012), molecules (Carr and Ye, 2009), ions (Blatt and Wineland, 2008; 
Schneider et al., 2012; Haffner et al., 2008) and photons (O'Brien, 2007), but also 
solid state systems (Clarke and Wilhelm, 2008; Wrachtrup and Jelezko, 2006; Han- 
son et al., 2007), we will summarize below recent advances in these directions with 
a particular focus on engineered open many-body systems and quantum simulations 
with atoms and ions. We note that in this review we intend, rather than providing 
a comprehensive overview of all recent developments in the field, to present our 
personal view on open-system quantum simulation, with a focus on work of the au- 
thors in Innsbruck in recent years ^ . Our emphasis is on presenting new concepts 
and building blocks, which we believe constitute first steps towards many-body sys- 
tems far away from thermodynamic equilibrium and future large-scale many-body 
simulations. 

Structure of this Review - Part 2 of this review presents theoretical and experimen- 
tal advances in digital quantum simulation with trapped ions and Rydberg atoms. In 
Sect. 2.1 we outline the basic concepts of digital quantum simulation and illustrate 
them by discussing results of recent experiments, which demonstrate the princi- 
ples of a digital quantum simulator in a trapped-ion quantum information processor 
(Lanyon et al., 2011) (Sect. 2.1). Subsequently, we discuss a proposal for a scal- 
able digital quantum simulator based on Rydberg atoms stored in optical lattices 
(Weimer et al., 2010). We show how this simulation architecture based on a multi- 
atom Rydberg gate (Miiller et al., 2009) allows one to simulate the Hamiltonian 
dynamics of spin models involving coherent «-body interactions such as Kitaev's 
toric code Hamiltonian (Sect. 2.2). In Sect. 2.3 we focus on digital simulation of 
dissipative many-body dynamics, which enables, e.g., the dissipative ground state 
preparation of the toric code via collective n-body dissipative processes. In this con- 
text, we discuss the corresponding reservoir-engineering techniques in the Rydberg 
simulator architecture, as well as recent experiments, which demonstrate the basic 
building blocks of an open-system quantum simulator with trapped ions (Barreiro 
et al., 201 1). Finally, we show how a combination of coherent and dissipative dy- 
namics might in the future enable the simulation of more complex spin models such 
as a three-dimensional U(l) lattice gauge theory. Finally, in Sect. 2.4 we comment 
on the effect of gate imperfections on the simulations. 

In part 3 we turn to engineered open many-body systems of cold atoms. In Sect. 3.1 
we demonstrate that quantum mechanical long-range order can be established dis- 
sipatively, and point out a route how this can be achieved via proper reservoir engi- 
neering, indeed extending the notion of quantum state engineering in cold atomic 



^ Parts of this review contain text and figure material from manuscripts by some of the 
authors, which have been published in other journals. 



gases from the Hamiltonian to the more general Liouvillian setting (Diehl et al., 
2008). We then give accounts for further central aspects of this general setting. 
In Sect. 3.2, we investigate the dynamical phase diagram resulting from the com- 
petition of unitary and dissipative dynamics, and identify several intrinsic many- 
body phenomena, underpinning that the stationary states of such systems consti- 
tute a novel class of artificial out-of-equilibrium ensembles (Diehl et al., 2010b). 
In Sect. 3.3, in the context of atomic fermions we reveal a novel dissipative pairing 
mechanism operative in the absence of any attractive forces (Diehl et al., 2010c), 
and point out how such systems may provide an attractive route towards quan- 
tum simulation of important condensed matter models, such as the Fermi-Hubbard 
model. Finally, we discuss in Sect. 3.4 how engineered dissipation may pave the 
way towards realizing in the lab topological states of matter (Diehl et al., 2011), 
and discuss some of their key many -body properties. 

We conclude with an outlook in Sect. 4, which summarizes present outstanding 
theoretical problems and challenges. 



2 Digital Quantum Simulation with Trapped Ions and Rydberg Atoms 



When is quantum simulation useful? As noted above, the main motivation for quan- 
tum simulation is to solve many-body problems where classical computers fail - 
or, at least, an efficient classical approach is presently not known. Indeed remark- 
able classical algorithms have been developed to solve specific problems and as- 
pects in equilibrium and out-of-equilibrium many body physics: examples include 
Monte-Carlo techniques (Ceperley, 1995; Prokof'ev et al., 1998), coupled-cluster 
expansion (Shavitt and Bartlett, 2009; Hammond et al., 1994), density functional 
theory (Parr and Yang, 1989), dynamical mean field theories (Georges et al., 1996), 
and density matrix renormalization group (DMRG) (SchoUwock, 2005; Hallberg, 
2006). These techniques may fail, when one encounters, for example, sign prob- 
lems in the Monte Carlo simulation of fermionic systems, or also in time depen- 
dent problems. An example is provided by quench dynamics: recent optical lattice 
experiments (Trotzky et al., 2011) have studied the time evolution after a quench, 
and a comparison with time-dependent DMRG calculations revealed the difficulty 
of predicting the long-time evolution due to growth of entanglement. These recent 
developments, enabled by the remarkable level of control achieved in analog cold- 
atom quantum simulators, are exciting, as they indicate for the first time possible 
large-scale entanglement in many-particle dynamics, close to the heart of quantum 
simulation. In the following section, we will outline the complementary route to 
simulate the time dynamics of interacting many-particle systems by the digital, i.e. 
gate-based, quantum simulation approach. 



2. 1 Concepts and First Experiments with Trapped Ions 



2.1.1 The Digital Simulation Method 



We start our discussion with the simulation of purely coherent dynamics gener- 
ated by a possibly time-dependent many-body Hamiltonian H{t) = Yicy Ha{t), and 
proceed in Sect. 2.3 with a detailed discussion of the digital simulation of dis- 
sipative dynamics according to many-body master Eqs. (2). It has been shown 
that a digital quantum simulator can implement the unitary time evolution op- 
erator U{t) generated by H{t) efficiently for any local quantum system (Lloyd, 
1996; Abrams and Lloyd, 1997; Ortiz et al., 2001), i.e., where the individual terms 
Ha are quasi-local. This means that they operate on a finite number of particles, 
due to interaction strengths that fall off with distance, for example. In this case 
it is possible to divide the simulation time t into small time steps A;^ = t/n and 
to implement the time evolution through a Trotter expansion of the propagator, 
U{t) - 0™=! ey^p{-iH(mAt)At). The key idea of the Trotter expansion is to approx- 
imate each propagator for a small time step according to the full Hamiltonian H{t) 
by a product of evolution operators for each quasi-local term, exp(-iH(mAt)At) ^ 
Y\a ex^p(-iHaimAt)At). In a digital quantum simulator each of the quasi-local prop- 
agators exp(-iHa(mAt)At) can be efficiently approximated by (or in many cases ex- 
actly decomposed into) a fixed number of operations from a universal set of gates 
(Lloyd, 1995; Kitaev, 1997; Nielsen and Chuang, 2000). As a consequence, the 
evolution is approximated by a stroboscopic sequence of many small time steps 
of dynamics due to the quasi-local interactions present in the system. The desired 
global time evolution according to the full many -body Hamiltonian,/? = -i[H{t),p] 
(see coherent part of Eq. (1)) emerges as an effective, coarse-grained description 
of the dynamics, as sketched in Fig. 1. For a finite number of time steps n, er- 
rors from possible non-commutativity of the quasi-local terms in the Hamiltonian, 
[Ha, Ha'] 4^ 0, are bounded (Nielsen and Chuang, 2000; Berry et al., 2007; Bravyi 
et al., 2008) and can be reduced by resorting to shorter time steps At or higher-order 
Trotter expansions (Suzuki, 1992). 



We note that it had been recognized early-on that dissipative dynamics can be ef- 
ficiently simulated by carrying out unitary dynamics on an enlarged Hilbert space 
(Lloyd, 1996), such that efficient simulation of Hamiltonian dynamics is in princi- 
ple sufficient to also realize open-system dynamics. In Sect. 2.3, we will discuss an 
alternative approach for digital simulation of dissipative dynamics, which combines 
unitary operations and dissipative elements (in a Markovian setting). Recently, ex- 
plicit error bounds for dissipative Trotter dynamics according to many-body master 
Eqs. (2) have been derived in Kliesch et al. (201 1). 
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Fig. 1. (Color online) Schematics of the working principle of a digital quantum simulator: 
For a specific (many-body) quantum system of interest to be simulated, the initial quan- 
tum state is stored in a register of qubits, which are encoded for instance in (meta-)stable 
electronic states of cold atoms in optical lattices or trapped ions. Then the time evolution 
of the system up to a time t is represented as a sequence of single- and many-qubit gates, 
according to a Trotter decomposition of the time evolution operator for small time steps 
At. Thus, the effective dynamics according to the desired model Hamiltonian //gff arises 
approximately and on a coarse-grained time scale. This digital, i.e., gate-based simulation 
approach is very flexible as the simulated (w-body) interactions can be substantially dif- 
ferent from and more complex than the physical one- and two-body interactions, which 
underlie the specific simulator architecture. The concept of digital quantum simulation is 
not limited to purely coherent Hamiltonian dynamics, but can be extended to the simulation 
of dissipative dynamics, as e.g. described by a many-body quantum master equation of the 
form of Eq. (1) with Liouvillian part of Eq. (2), and as discussed in detail in Sect. 2.3. 

2.1.2 Coherent Digital Simulation with Trapped Ions 

A recent experiment carried out on a small-scale trapped ion quantum computer 
(Lanyon et al., 2011) has explored and demonstrated in the laboratory the various 
aspects of digital Hamiltonian quantum simulation. In a series of digital quantum 
simulations according to different types of interacting quantum spin models the 
performance of the digital simulation approach for systems of increasing complex- 
ity in the interactions and increasing system sizes was quantitatively studied. The 
experiments, whose main aspects we will briefly summarize in this section, have 
been enabled by remarkable progress in the implementation of individual gate oper- 
ations (see Fig. 2 for details on the experimental simulation toolbox). In particular, 
multi-ion entangling gates have been realized with fidelities higher than 99% for 
two ions (Benhelm et al., 2008; Roos, 2008), and for up to 14 qubits (Monz et al., 
2011). 

Exploring Trotter dynamics with two spins - To illustrate the Trotter simulation 
method, the conceptually most simple example of an Ising system of two inter- 
acting spin- 1/2 particles as an elementary building block of larger and more com- 
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Fig. 2. (Color online) Toolbox of quantum operations in the Innsbruck ion trap quantum 
computer, a) Simplified level scheme of laser-cooled "^''Ca^ ions stored in a linear Paul 
trap: Long-lived internal electronic states \D} = |0> and \S) = |1> represent the qubit, while 
short-hved transitions are used for read-out of the quantum state of the qubit using a flu- 
orescence measurement technique, b) The universal set of gates is formed by addressed 
single-qubit z-rotations and c) collective x- and y-rotations as well collective entangling 
operations U^i , as suggested by M0lmer and S0rensen (1999). d) For the simulation of 
open-system dynamics (see Sect. 2.3) the string of ions can be divided into system qubits 
S (ions 1 through n) and an "environment" qubit E. Coherent gate operations on S and E, 
combined with a controllable dissipative mechanism involving spontaneous emission of a 
photon from the environment ion via an addressed optical pumping technique (Schindler 
et al., 201 1), allow one to tailor the coupling of the system qubits to an artificial environ- 
ment (see Barreiro et al. (2011) for experimental details). This should be contrasted to 
the residual, detrimental coupling of the system (and environment) ions to their physical 
environment. Figure adapted from Barreiro et al. (2011). 

plex spin models was studied: The Hamiltonian is given by the sum of two non- 
commuting terms, H = Hint + ^magn^ where H^nt = Jc\crl describew a spin- 
spin interaction, and //magn the coupling to an effective, transverse magnetic field 
^magn = ^(cj + o"^)- This was one of the first systems to be simulated with trapped 
ions following an analog approach (Friedenauer et al., 2008; Porras and Cirac, 
2004). The experiments (Lanyon et al., 2011) studied the two-spin dynamics both 
for the time-independent Ising-Hamiltonian (see Fig. 3(a)), as well as for the time- 
dependent case where the interaction term Hi^t was linearly ramped up in time (see 
Fig. 3(b)). The time evolution was realized by a first-order Trotter decomposition. 
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Fig. 3. (Color online) Digital Hamiltonian simulation with trapped ions. This figure and 
the following one present some basic concepts of digital Hamiltonian simulation, and illus- 
trate them with examples from a recent experiment with trapped ions (Lanyon et al., 201 1), 
where the digital approach was used to simulate various interacting quantum spin models 
of different complexity in the interactions and different system sizes. The simulations were 
realized using the toolbox of available coherent gates specified in Fig. 2. a) Time-indepen- 
dent Hamiltonian simulation. Dynamics of the initial state |tT) under a time-independent 
two-spin Ising Hamiltonian with J = 2B: As expected, the simulated dynamics according 
to a first-order Trotter decomposition converge closer to the exact dynamics as the digital 
resolution is increased, i.e. the size of the individual time steps is decreased. It is conve- 



nient to introduce a dimensionless Hamiltonian H, i.e. H=EH such that U=e 



_„-iHEM 



and 



the evolution is quantified by a unitless phase = EAt. Each single digital step is given 
by U1U2 - U^2{6a/n) U (i.2){9aln) with 9a = n/lyl and n - \ and n = 4 (finer Trotter 
resolution). (Labeling: Lines: exact dynamics. Unfilled shapes: ideal digitised (Trotter de- 
composition). Filled shapes: experimental data. Bft, ♦J.J-)- b) Time -dependent Hamiltonian 
dynamics. Time evolution under a two-spin Ising Hamiltonian, where the spin-spin inter- 
action strength J increases linearly from to 4B during a total evolution given by 9t=7i/2. 
In the experiment, the continuous dynamics is approximated using a sequence of 24 gates, 
with c=U (i.2)(7r/8), d=Ug2{n/l6). The increase of / over time is reflected by an increase in 
the number of J-blocks per Trotter step. The observed oscillation in population expectation 
values (measured in the cTji-basis) is a diabatic effect due to the finite speed in ramping up 
the interaction term Him (•^^x. A*— <— x)- Percentages: fidelities between measured and 
exact states with uncertainties less than 2%. Figure reprinted with permission from Lanyon 
et al. (201 1). Copyright 201 1 by MacMillan. 
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Fig. 4. (Color online) Continuation of digital Hamiltonian simulation with trapped ions. 
Simulation ofn-body interactions, a) In a digital simulation, «-body spin interactions (with 
n > 2) are usually realized by quantum circuits involving 2m two-qubit C-NOT gates 
(Nielsen and Chuang, 2000). However, the availability of high-fidelity, collective entan- 
gling gates acting on n ions allows one to bundle the effect of such series of two-qubit 
gates and thus to realize, e.g., six-body interactions by a highly compact, experimentally 
efficient quantum circuit involving two six-ion gates {4D=Ug2{7T/4)), interspersed with one 
single-ion z-rotation {F=U (i)(20)). The strength of the six-body interaction is controlled 
by the phase 6 in the single-qubit rotation (see Ref. Miiller et al. (201 1) for theoretical de- 
tails), b) Experimentally observed dynamics induced by a six-body spin interaction, which 
directly couples the states ItTTTTT) and UJ,J,J.J,J,), periodically producing a maximally en- 
tangled GHZ state. Lines: exact dynamics. Filled shapes: experimental data. The quanti- 
tative characterization and assessment of errors of such multi-qubit building blocks is a 
non-trivial task, as standard quantum process tomography is impractical for more than 3 
qubits. The inequality bounds the quantum process fidelity Fp at 9=0.25 - see online ma- 
terial of Lanyon et al. (2011) for details on the employed technique. Figure reprinted with 
permission from Lanyon et al. (20 11). Copyright 20 11 by MacMillan. 

where the propagators for small time steps according to the two Hamiltonian terms 
were decomposed into sequences of experimentally available single- and two-qubit 
gates. 

Simulation of larger systems and n-body interactions - Experiments with up to six 
ions (Lanyon et al., 2011) showed that the digital approach allows arbitrary in- 
teraction distributions for larger interacting spin systems to be programmed. For 
instance, it is possible to implement spatially inhomogeneous distributions of in- 
teraction strengths and to simulate n-body interaction terms, with n > 2, in a 
non-perturbative way (see Fig. 4). Many-body spin interactions of this kind are 
an important ingredient in the simulation of systems with strict symmetry require- 
ments (Kassal et al., 2011). Furthermore, they appear in the context of many-body 
quantum systems exhibiting topological order (Nayak et al., 2008) and in the con- 
text of topological quantum computing and memories (Kitaev, 2003; Dennis et al., 
2002). In Sect. 2.2.1 we will discuss in more detail Kitaev's toric model (Kitaev, 
2003) as an example for a complex spin model involving four-body spin interac- 
tion terms. Engineering of three-body interactions in analog quantum simulators 



11 



has been suggested for trapped ions (Bermudez et al., 2009) and polar molecules 
(Biichler et al., 2007); however, it is in general very difficult to achieve dominant, 
higher-order interactions of substantial strength via analog quantum simulation 
techniques. Fig. 4(b) shows the digital simulation of time evolution according to 
a six-spin many-body interaction, where each Trotter time step was experimentally 
realized by a highly compact quantum circuit involving two collective six-ion en- 
tangling gates as essential resource (Lanyon et al., 201 1). 

In view of these remarkable experimental advances and the demonstrated flexibil- 
ity and control achieved so far, two major remaining challenges are (i) the quan- 
tum simulation of open-system quantum dynamics according to many-body master 
equations of the form (1) and (ii) to scale up the simulations from a few qubits to 
larger system sizes. Regarding the latter aspect, we will in the next section leave 
the trapped ions for a moment and switch to another physical platform, where we 
will discuss an a priori scalable, digital simulation architecture based on Rydberg 
atoms stored in optical lattices or magnetic trap arrays. In Sect. 2.3 we will then 
extend the discussion to open many-particle quantum systems and describe how to 
simulate complex dissipative many -body dynamics. In this context we will come 
back to trapped ions, where recently the building blocks of an open-system quan- 
tum simulator have been successfully implemented (Barreiro et al., 201 1). 



2.2 Scalable Quantum Simulation with Rydberg Atoms 



Laser excited Rydberg atoms (Gallagher, 1994) offer unique possibilities for quan- 
tum information processing and the study of strongly correlated many-body dynam- 
ics. Atoms excited to high-lying Rydberg states interact via strong and long-range 
dipole-dipole or Van der Waals forces (Gallagher, 1994) over distances of several 
yum, which are internal state-dependent and can be up to 12 orders of magnitude 
stronger than interactions between ground state atoms at a comparable distance 
(Saffman et al., 2010). Electronic level shifts associated with these interactions can 
be used to block transitions of more than one Rydberg excitation in mesoscopic 
atomic ensembles. This "dipole blockade" (Jaksch et al., 2000; Lukin et al., 2001) 
mechanism underlies the formation of "superatoms" in atomic gases with a sin- 
gle Rydberg excitation shared by many atoms within a blockade radius. This ef- 
fect gives rise to strongly correlated, dominantly coherent many-body dynamics 
(Raitzsch et al., 2008), which has been explored in recent years both experimen- 
tally (Tong et al., 2004; Singer et al., 2004; Cubel et al., 2005; Vogt et al., 2006; 
Mohapatra et al., 2007; Heidemann et al., 2007; Reetz-Lamour et al., 2008) and 
theoretically (Pohl et al., 2010; Weimer et al., 2008; Olmos et al., 2009; Sun and 
Robicheaux, 2008; Honer et al., 2010). In the context of quantum information pro- 
cessing, it has been recognized that these strong, switchable interactions between 
pairs of atoms potentially provides fast and addressable two-qubit entangling op- 
erations (Jaksch et al., 2000; Lukin et al., 2001) or effective spin-spin interactions 
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(Lesanovsky, 2011; Pohl et al., 2010; Schachenmayer et al., 2010; Weimer et al., 
2008); recent theoretical proposals have extended Rydberg -based protocols towards 
a single-step, high-fidelity entanglement of a mesoscopic number of atoms (M0ller 
et al., 2008; Miiller et al., 2009). Remarkably, the basic building blocks of Rydberg- 
based quantum information processing have been demonstrated recently in the lab- 
oratory by several groups, which observed the dipole blockade between a pair of 
neutral Rydberg atoms stored in optical tweezers (Urban et al., 2009; Gaetan et al., 
2009). Here, the Rydberg blockade was used as a mechanism to create two-atom 
entanglement (Wilk et al., 2010) and to realize the first neutral atom two-qubit C- 
NOT gate (Isenhower et al, 2010). 

On the other hand, cold atoms stored in optical lattices or magnetic trap arrays of- 
fer a versatile platform for a priori scalable quantum information processing and 
quantum simulation (Jaksch et al., 1998a; Jaksch and ZoUer, 2005; Greiner et al., 
2001; Bloch et al., 2008; Lewenstein et al., 2006; Dalibard et al., 2011). In partic- 
ular, in sufficiently deep lattices, where tunneling between neighboring lattice sites 
is suppressed, single atoms can be loaded and kept effectively frozen at each lat- 
tice site, with long-lived atomic ground states representing qubits or effective spin 
degrees of freedom. Working with large-spacing lattices, with inter-site distances 
of the order of a few //m (Nelson et al., 2007; Whitlock et al., 2009) allows single- 
site addressing with laser light, and thus individual manipulation and readout of 
atomic spins. Very recently, several groups have achieved single-site addressing in 
optical lattices (Bakr et al., 2010a; Sherson et al., 2010; Bakr et al., 2010b) and 
manipulation of individual spins in this setup (Weitenberg et al., 201 1) (see Fig. 5). 



As we will discuss below, given these achievements and the future integration of 
techniques for coherent laser excitation of Rydberg atoms in addressable (optical) 
lattice setups (Viteau et al., 2011; Anderson et al., 2011), in principle all essen- 
tial ingredients seem to exist already in the laboratory to build a scalable, digital 
quantum simulator based on cold Rydberg atoms (Weimer et al., 2010). 



Before specifying in more detail the concrete physical architecture of the Rydberg 
quantum simulator proposed in Weimer et al. (2010), we will in the next section 
discuss a specific many-body spin model of interest: Kitaev's toric code (Kitaev, 
2003). This model represents a paradigmatic example of a large class of spin mod- 
els, which have in the last years attracted great interest in the context of topological 
quantum information processing and as strongly interacting many-body quantum 
systems exhibiting topological order (Nayak et al., 2008; Wen, 2004). This exam- 
ple illustrates the generic challenges and goals of a quantum simulation of complex 
many-body models, which are to be addressed in a concrete physical implementa- 
tion of a quantum simulator. The realization of a more complex setup of a three- 
dimensional t/(l) lattice gauge theory giving rise to a spin liquid phase will be 
discussed below in Sect. 2.3.3. 
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Addressing laser beam 





Atoms in two-dimensional optical iattice 




Fig. 5. (Color online) Single-site addressing of atoms in an optical lattice. The left part 
shows a schematics of atoms loaded into a square optical lattice, where they form a Mott 
insulator state with one atom per lattice site. Atoms residing on individual lattice sites in the 
x-y plane can be optically addressed with an off -resonant laser beam, which can be focused 
to individual sites by means of a high-aperture microscope objective. The upper part of (a) 
shows an experimentally obtained fluorescence image of a Mott insulator site with one 
atom per site, where a subset of atoms (diagonal of the image) has been transferred from 
an internal state |0) to 1 1 ) by means of the single-site addressed beam. Before fluorescence 
detection, the atoms in |1) are removed from the lattice by a resonant laser pulse. The 
bottom part shows the reconstructed atom number distribution (see Sherson et al. (2010) 
for details on the reconstruction algorithm), where filled black circles correspond to single 
atoms and dots indicate the position of the lattice sites. Figure adapted with permission 
from Weitenberg et al. (2011). Copyright 2011 by MacMillan. 

2.2.1 Paradigmatic Example: Simulation ofKitaev's toric Code Hamiltonian 

Kitaev's toric code is a paradigmatic, exactly solvable model, out of a large class of 
spin models, which have recently attracted a lot of interest in the context of stud- 
ies on topological order and quantum computation. It considers a two-dimensional 
setup, where spins are located on the edges of a square lattice (Kitaev, 2003). The 
Hamiltonian H = -Eq (Z^^p + Z.« BA is a sum of mutually commuting stabilizer 
operators Ap = Yliep crj and Bg = Ylies 0"^ which describe four-body interactions 
between spins located around plaquettes (Ap) and vertices (5,) of the square lattice 
(see Fig. 6a). All Ap and B^ stabilizer operators mutually commute, thus the ground 
state of the Hamiltonian is a simultaneous eigenstate of all stabilizer operators Ap 
and Bs with eigenvalues -i-l, and gives rise to a topological phase: the ground state 
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Fig. 6. (Color online) Schematics of the Rydberg quantum simulator architecture and a 
multi-atom C-NOT^ Rydberg gate as its principal building block, a) The Rydberg quantum 
simulator (Weimer et al., 2010) is particularly suited for the simulation of coherent and dis- 
sipative dynamics of complex quantum spin models involving «-body interactions and con- 
straints. A paradigmatic example is Kitaev's toric code Hamiltonian (Kitaev, 2003), where 
spins are located on the edges of a two-dimensional square lattice and interact via four- 
body plaquette or vertex interactions. The model exhibits two types of localized quasi-par- 
ticle excitations (depicted as red and green dots), which exhibit Abelian anyonic statistics 
under braiding, i.e. when they are winded around each other, b) A mesoscopic multi-atom 
Rydberg gate (Miiller et al., 2009) applied to subsets of four spins around plaquettes and 
vertices, and additional control atoms, which are located at the centers of the plaquettes and 
on the vertices of the lattice, allows one to efficiently realize such many-body plaquette and 
vertex interactions. Here, controllable strong and long-range Rydberg interactions mediate 
effective four-body interactions among the system spins. By a combination of the multi- 
-qubit C-NOT gate shown in (c) with optical pumping on the auxiliary control atoms, it is 
possible to engineer dissipative «-body processes. This many-body reservoir engineering 
can be used to realize cooling dynamics, which leads, e.g., to the dissipative ground state 
preparation of Kitaev's toric code Hamiltonian. 

degeneracy depends on the boundary conditions and topology of the setup, and the 
elementary excitations exhibit Abelian anyonic statistics under braiding, i.e. when 
they are winded around each other. The toric code shows two types of localized ex- 
citations corresponding to -1 eigenstates of each stabilizer Ap ("magnetic charge", 
filled red dots in Fig. 6a) and Bp ("electric charge", filled green dots). 

In addition to the toric code Hamiltonian, one can formulate a dissipative many- 
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body dynamics, which "cools" into the ground state manifold of the many-body 
Hamiltonian. Such dissipative time evolution is provided by a Liouvillian (2) with 
quantum jump operators, 

Cp = -(T]{\-Ap), c, = -(r){,\-B,), (3) 

with / 6 p and j 6 s, which act on four spins located around plaquettes p and 
vertices s, respectively. In Sect. 2.3.2 we will discuss in detail how these four-body 
quantum jump operators can be physically implemented in the Rydberg simulator 
architecture of Weimer et al. (2010). The jump operators are readily understood 
as operators which "pump" from -1 into -i-l eigenstates of the stabilizer operators: 
the part (1 - Ap)/2 of Cp is a projector onto the eigenspace of A^ with -1 eigen- 
value (an excited state with a "magnetic charge" present), while all states in the -l-l 
eigenspace are dark states. The subsequent spin flip 07 transfers the excitation to 
the neighboring plaquette. The jump operators then give rise to a random walk of 
anyonic excitations on the lattice, and whenever two excitations of the same type 
meet they are annihilated, resulting in a cooling process, see Fig. 7a. Similar argu- 
ments apply to the jump operators Cs- Efficient cooling is achieved by alternating 
the index i of the spin, which is flipped. 

Our choice of the jump operator follows the idea of reservoir engineering of inter- 
acting many-body systems as discussed in Diehl et al. (2008); Kraus et al. (2008) 
and in Sect. 3. In contrast to alternative schemes for measurement based state prepa- 
ration (Aguado et al., 2008), here, the cooling is part of the time evolution of the 
system. These ideas can be readily generalized to more complex stabilizer states 
and to setups in higher dimensions, as in, e.g., the color codes developed in Bombin 
and Martin-Delgado (2006, 2007), and the simulation of a three-dimensional U(l) 
lattice gauge theory, which will be discussed in Sect. 2.3.3. 

In conclusion, the main challenge in the quantum simulation of coherent Hamilto- 
nian dynamics and dissipative ground state preparation of many-body spin models 
such as Kitaev's toric code Hamiltonian lies in (i) the realization of strong n-body 
interactions, and (ii) the ability to tailor multi-particle couplings of the many-body 
system to a reservoir, such that the dissipative dynamics gives rise to ground state 
cooling, as described by a many-body master Eq. (2) with many-body quantum 
jump operators of Eq.(3). 



2.2.2 A Mesoscopic Rydberg Gate 

Let us now turn to the physical implementation of the digital Rydberg simulator 
setup suggested in Weimer et al. (2010). A key ingredient of the proposed archi- 
tecture are additional auxiliary qubit atoms in the lattice, which play a two-fold 
role: First, they control and mediate eff'ective n-body spin interactions among a 
subset of n system spins residing in their neighborhood in the lattice, as e.g. the 
four-body plaquette and vertex interactions of Kitaev's toric code Hamiltonian dis- 
cussed above. In the proposed scheme this is achieved efficiently making use of 
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Fig. 7. (Color online) Cooling of Kitaev's toric code: a) A dissipative time step incoherently 
moves one anyonic excitation (red dot) on top of a second anyon located on a neighboring 
plaquette, annihilating each other and thus lowering the internal energy of the system. The 
anyon of the other type (an "electric charge", filled green dot located on a vertex of the 
lattice) remains unaffected by this cooling step, b) Numerical simulation of the cooling for 
N lattice sites (periodic boundary conditions). Single trajectories for the anyon density n 
over time are shown as solid lines. Filled circles represent averages over 1000 trajectories. 
The initial state for the simulations is the fully polarized, experimentally easily accessible 
state of all spins down. For perfect gates in the digital quantum simulation discussed in 
detail in Sect. 2.3.2, the energy of the system reaches the ground state energy in the long 
time limit, while for imperfect gates heating events can occur (blue solid line) and a finite 
density of anyons n remains present (blue circles). The characteristic time scale k~^ for 
cooling is set by (i) the gate parameters in the quantum circuit decomposition discussed 
below (see Sect. 2.3.2 and (ii) by the duration for the implementation of the underlying 
quantum gates. Figure reprinted with permission from Weimer et al. (2010). Copyright 
2010 by MacMillan. 

single-site addressability and a parallelized multi-qubit gate, which is based on a 
combination of strong and long-range Rydberg interactions and electromagneti- 
cally induced transparency (EIT) and is schematically shown in Fig. 6b. This gate 
has been suggested and analyzed in Miiller et al. (2009). As it plays a central role 
in the simulation architecture, we will briefly and on a qualitative level review its 
main features here. Second, the auxiliary atoms can be optically pumped, thereby 
providing a dissipative element, which in combination with Rydberg interactions 
results in effective collective dissipative dynamics of a set of spins located in the 
vicinity of the auxiliary particle. This enables, e.g., the simulation of dissipative 
dynamics for ground state cooling of Kitaev's toric code and related models. 

Setup of the Rydberg gate - The envisioned setup is illustrated in Fig. 6b. A control 
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atom and a mesoscopic ensemble of, say, four atoms are stored in separate trapping 
potentials, e.g. in two dipole traps as in Wilk et al. (2010); Isenhower et al. (2010) 
or in neighboring lattice sites of a (large-spacing) optical lattices or magnetic trap 
array (Whitlock et al., 2009). The multi-qubit gate exploits state-dependent Ryd- 
berg interactions and reaUzes a controUed-NOT'^ (CNOT^) gate, which is defined 
by 

[/, = |0)<0|,Sl, + |l)<l|,S<. (4) 

('= 1 (■= 1 

Depending on the state of the control qubit - the state of all N target qubits is left 
unchanged or flipped. Here, |0), |1) and |A), \B) denote long-lived ground states of 
the control and ensemble atoms, respectively (see Fig. 6b), and cr^\A)i = |5), and 

The basic elements of the gate of Eq. (4) are: (i) the control atom can be individu- 
ally addressed and laser excited to a Rydberg state conditional to its internal state, 
thus (ii) turning on or off" the strong long-range Rydberg-Rydberg interactions of 
the control with ensemble atoms, which (iii) via EIT-type interference suppresses 
or allows the transfer of all ensemble atoms from \A) or \B) conditional to the state 
of the control atom. It does not necessarily require individual addressing of the en- 
semble atoms, in contrast to a possible implementation of the gate (4) by a sequence 
of A^ two-qubit C-NOT gates. 

Implementation of the gate operation - For the physical realization of the operation 
(4), an auxiliary Rydberg level \r) of the control atom is used, which is resonantly 
coupled to 1 1) by a laser with (two-photon) Rabi frequency Qr (see Fig. 8). For 
the ensemble atoms the two stable ground states \A) and \B) are coupled far off"- 
resonantly in a A-configuration with Rabi frequency Qp and detuning A to a low- 
lying, intermediate state \P) (e.g. 5^/'3/2 in case of *^^Rb). A second laser with Rabi 
frequency Q^ (A » Q^ > ^p) couples \P) to a Rydberg state \R) of the ensemble 
atoms, such that the two ground states |A) and \B) are in two-photon resonance with 
|i?), as depicted in Fig. 8b. 

The conditional, coherent transfer of population between the ground states of the 
ensemble atoms, as required for the C-NOT operation, is then achieved by a se- 
quence of three laser pulses (shown in Fig. 8a): (i) a short 7r-pulse on the control 
atom, (ii) a smooth Raman 7r-pulse ^pit) with ^ dtQp(t)/(2A) = n acting on all 
ensemble atoms, and (iii) finally a second ;r-pulse on the control atom. The effect 
of this pulse sequence can be understood by distinguishing the two possible cases 
of (a) blocked transfer (for the control atom initially residing in the logical state |0)) 
and (b) enabled transfer (with the control atom initially in |1)). 

(a) Blocked population transfer: For the control atom initially residing in the logical 
state |0) the first pulse has no effect. In the regime Qp «; Q.c, the laser configuration 
of the ensemble atoms realizes an FIT scenario (Fleischhauer et al., 2005), where 
the strong always-on "control" laser field Q.^ suppresses via destructive interfer- 
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Fig. 8. (Color online) Mesoscopic Rydberg gate, a) Sequence of laser pulses (not to scale), 
b) Electronic level structure of the control and ensemble atoms. The ground state |1) is 
resonantly coupled to the Rydberg state \r). The states |A) and |B) are off-resonantly coupled 
(detuning A, Rabi frequency Op) to \P). A strong laser with Rabi frequency Qc ^ ^p 
couples the Rydberg level \R) to \P) such that \R) is in two-photon resonance with \A) and 
\B). In this situation (known as EIT) Raman transfer from \A) to \B) is inhibited, c) With 
the control atom excited to \r) the two-photon resonance condition is lifted as the level \R) 
is shifted due to the interaction energy V between the Rydberg states, thereby enabling 
off-resonant Raman transfer from \A) to \B). Figure adapted from Miiller et al. (2009) 



ence coupling of the "probe" laser Qp to the intermediate state \P) and thus also the 
second-order Raman coupling. This also effectively inhibits population transfer be- 
tween \A) and \B). As a consequence, the Raman pulses Qp are ineffective (as well 
as the second 7r-pulse applied to the control atom in |1)). The state of the ensemble 
atoms remains unchanged, thereby realizing the first logical half of the gate (4). 

(b) Enabled population transfer: If the control atom initially resides in |1), it is 
excited to the Rydberg state \r) by the first pulse. Due to strong repulsive Rydberg 
interactions V > between the control atom in \r) and ensemble atoms in \R), the 
Rydberg level of the ensemble atom is now shifted by the energy V (see Fig. 8c), de- 
spite the fact that the Rydberg state \R) of the ensemble atoms is not populated. This 
interaction-induced energy shift lifts the two-photon resonance condition, which 
underlies the EIT scenario and is crucial to block the Raman transfer between \A) 
and \B). Now, the Raman lasers couple ofF-resonantly to IP) and the coherent pop- 
ulation transfer between |A) and \B) takes place. 
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A quantitative analysis of the gate performance in Miiller et al. (2009) shows that 
the effect of the relevant error sources such as radiative decay from the \P) and the 
Rydberg states and possible mechanical effects are negligible for realistic atomic 
and laser parameters. Remarkably, undesired destructive many-body effects orig- 
inating from undesired, but possibly strong Rydberg interactions between the en- 
semble atoms can be effectively suppressed and minimized in the limit Qp «c Qc- 
As a consequence, the gate also works reliably and with high fidelity for a moder- 
ate number of ensemble atoms separated by up to a few microns, it is robust with 
respect to inhomogeneous inter-particle distances and varying interaction strengths 
and can be carried out on a microsecond timescale (Miiller et al., 2009). 



2.2.3 Simulation of Coherent Many -Body Interactions 

The many-qubit Rydberg gate (Miiller et al., 2009) discussed in the previous section 
is the key building block of the Rydberg quantum simulator architecture (Weimer 
et al., 2010). Using an auxiliary qubit located at the center of a four-atom pla- 
quette allows one to efficiently simulate coherent n-body interactions such as the 
four-body spin plaquette interactions Ap = O; cr^ appearing in Kitaev's toric code 
Hamiltonian (Fig. 9). The general approach consists of a sequence of three coher- 
ent steps, as depicted in Fig. 9b: (i) First, a gate sequence M is performed, which 
coherently encodes the information whether the four system spins are in a -l-l or -1 
eigenstate of Ap in the two internal states of the auxiliary control qubit (see Fig. 9c). 
(ii) In a second step, a single qubit-gate operation, which depends on the internal 
state of the control qubit, is applied. Due to the previous mapping this manipulation 
of the control qubit is equivalent to manipulating the subspaces with fixed eigenval- 
ues ±1 of Ap. Thus, effectively, the application of a single-qubit gate exp(-?'0cr^) 
on the control qubit imprints a phase shift exp(+?'0) on all ±1 eigenstates of the sta- 
bilizer operator Ap. (iii) Finally, the mapping M is reversed, and the control qubit 
returns to its initial state |0). Consequently, at the end of the sequence, the auxiliary 
qubit effectively factors out from the dynamics of the four system spins, which in 
turn have evolved according to the desired time evolution 

U = exp(-i(f)Ap) = M'^ exp (-/0cr^) M. (5) 

Note that the essential resource for one time step consists of two applications of the 
mesoscopic Rydberg gate U^, which up to local rotations realizes the mappings M 
and M~^. In contrast, a standard implementation via two-qubit C-NOT gates would 
correspond to eight entangling operations (Nielsen and Chuang, 2000). 

For small phase imprints «c 1 the mapping reduces to the standard equation for 
coherent time evolution according to the master equation dtP = -iEo[Ap,p] + oicfP') 
and thus implements the propagator for a small Trotter time step according to the 
four-body spin interaction Ap on one plaquette. The above scheme for the imple- 
mentation of the many-body interaction Ap can be naturally extended to arbitrary 
many-body interactions between the system spins surrounding the control atom, as 
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Fig. 9. (Color online) Simulation of coherent «-body interactions, a) Kitaev plaquette term 
corresponding to four-body interactions Ap - Y\t=i ^^- b) Three-step gate sequence, which 
implements desired time evolution U = exp(-/(/>Ap) of the four system spins, mediated 
by an auxiliary control qubit. c) The gate sequence M coherently maps the information, 

whether the system spins reside in a -i-l (e.g. \+ -\ » or -1 eigenstate (e.g. |-i i--i-)) of 

the many -body interaction Ap onto the internal state |0>c and 1 1 )c of the control qubit. The 
mapping is given by M = exp{+incrl/4) Ug exp(-/;r(7^/4), i.e., up to single-qubit y-ro- 
tations of the control qubit, by the multi-atom Rydberg gate of Eq. (4). After the map- 
ping, a single-qubit z-rotation of the control qubit exp {-icpcrl) effectively imprints a phase 
exp(+/0) on all ±1 eigenstates of Ap. After the mapping M is reversed, the control qubit re- 
turns to |0>c and thus factors out from the dynamics of the system spins, which have evolved 
according to U. 

e.g., the Bp interaction terms in the toric code. Gate operations on single system 
spins allow to transform a^. into (P. and 07, in accordance with previous proposals 
for digital simulation of spin Hamiltonians (S0rensen and M0lmer, 1999), while 
selecting only certain spins to participate in the many-body gate via local address- 
ability gives rise to the identity operator for the non-participating spins. 

The associated energy scale of the many-body interactions becomes Eq = 0/t with 
T the physical time needed for the implementation of all gates, which are required 
for a single time step according to the many-body Hamiltonian on the whole lat- 
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tice. Note that in principle many of these operations at sufficiently distant areas of 
the lattice can be done in parallel, for instance by using super-lattices (Lee et al., 
2007; Foiling et al., 2007) for the application of the required laser pulses. In this 
case the energy scale Eq becomes independent of the lattice size, and is essentially 
only limited by the fast micro-second time scale of the Rydberg gates, potentially 
allowing for characteristic energy scales Eq on the order of 10-100 kHz (Weimer 
etal., 2010, 2011). 



2.3 Digital Simulation of Open-System Dynamics 



In the previous sections, we have focused on the principles and physical exam- 
ples of digital simulation of coherent many-body interactions. Let us now extend 
the discussion to the digital simulation of dissipative many-body dynamics. The 
dynamics of an open quantum system S coupled to an environment E can be de- 
scribed by the unitary transformation psE ^-^ UpseU\ with psE the joint density 
matrix of the composite system S -I- E. Thus, the reduced density operator of the 
system will evolve as p = Txe{UPseU^)- The time evolution of the system can also 
be described by a completely positive Kraus map (Nielsen and Chuang, 2000) 

p ^ &{p) = ^ EkpEl, (6) 

k 

where p denotes the reduced density operator of the system, {E/,} is a set of op- 
eration elements satisfying YjkElEi, = 1, and we assume an initially uncorrelated 
system and environment. For the case of a closed system, decoupled from the en- 
vironment, the map of Eq. (6) reduces to p i-^ t/pt/^ with U the unitary time 
evolution operator of the system. The Markovian limit of the general quantum op- 
eration (6) for the coherent and dissipative dynamics of a many-particle system is 
given by the many-body master Eq. (1) discussed above. 

Control of both coherent and dissipative dynamics is then achieved by finding cor- 
responding sequences of maps specified by sets of operation elements {Et} and 
engineering these sequences in the laboratory. In particular, for the example of dis- 
sipative quantum state preparation, pumping to an entangled state |i/f) reduces to 
implementing appropriate sequences of dissipative maps. These maps are chosen 
to drive the system to the desired target state irrespective of its initial state. The re- 
sulting dynamics have then the pure state |i/r) as the unique attractor, p i-^ |iA){iAl- In 
quantum optics and atomic physics, techniques of optical pumping and laser cool- 
ing are successfully used for the dissipative preparation of quantum states, although 
on a single-particle level. The engineering of dissipative maps for the preparation 
of entangled states can be seen as a generalization of this concept of pumping and 
cooling in driven dissipative systems to a many-particle context. For a discussion of 
Kraus map engineering from a control-theoretical viewpoint see also the literature 
(Lloyd and Viola, 2001; Wu et al., 2007; Bolognani and Ticozzi, 2010; Verstraete 
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et al., 2009) and the discussion on open - vs. closed-loop simulation scenarios at 
the end of Sect. 2.3.2. To be concrete, here we focus on dissipative preparation of 
stabilizer states, which represent a large family of entangled states, including graph 
states and error-correcting codes (Steane, 1996; Calderbank and Shor, 1996). Sim- 
ilar ideas for dissipative preparation of correlated quantum phases are discussed in 
Sect. 3 in the context of analog many-body quantum simulation in cold bosonic 
and fermionic atomic systems. 



2.3.1 Bell State Pumping 

Before discussing the dissipative preparation of many-body phases such as ground 
state cooling of Kitaev's toric code Hamiltonian, we start by outlining the con- 
cept of dissipative Kraus map engineering for the simplest non-trivial example 
of "cooling" a system of two qubits into a Bell state. The Hilbert space of two 
qubits is spanned by the four Bell states defined as |0*) = -^(|00) + |11)) and 

|*P*) = -4=(|01) + |10)). Here, |0) and |1) denote the computational basis of each 
qubit, and we use the short-hand notation |00) = |0)i|0)2, for example. These max- 
imally entangled states are stabilizer states: the Bell state |0^), for instance, is said 
to be stabilized by the two stabilizer operators Z1Z2 and X1X2, where X and Z de- 
note the usual Pauli matrices, as it is the only two-qubit state being an eigenstate 
of eigenvalue -i-l of these two commuting observables, i.e. ZiZ2|<l)^) = |0^) and 
X1X2IO+) = |0^). In fact, each of the four Bell states is uniquely determined as 
an eigenstate with eigenvalues ±1 with respect to Z1Z2 and X1X2. The key idea of 
cooling is that we can achieve dissipative dynamics which pump the system into 
a particular Bell state, for example p \-^ |*P")(*P"|, by constructing two dissipative 
maps, under which the two qubits are irreversibly transfered from the -l-l into the 
-1 eigenspaces of Z1Z2 and X1X2, as sketched in the upper part of Fig. 10. The dis- 
sipative maps are engineered with the aid of an ancilla "environment" qubit (Lloyd 
and Viola, 2001 ; Diir et al., 2008) and a quantum circuit of coherent and dissipative 
operations. 

Kraus maps for Bell state pumping - For Z1Z2, the dissipative map which induces 
pumping into the -1 eigenspace is given by p h^ S(p) = EipE\ + E2pE\ with 

Ei= Vp^2^(1+ZiZ2), 

£2 = ^(l-ZiZ2)+ Vi^^(l+ZlZ2). (7) 

The map's action as a uni-directional pumping process can be seen as follows: 
since the operation element Ei contains the projector ^(1 -I- Z1Z2) onto the -i-l 
eigenspace of Z1Z2, the spin flip X2 can then convert -l-l into -1 eigenstates of Z1Z2, 
e.g., |0+) 1-^ \^^). In contrast, the -1 eigenspace of Z1Z2 is left invariant. The cool- 
ing dynamics are determined by the probability of pumping from the -i-l into the -1 
stabilizer eigenspaces, which can be directly controlled by varying the parameters 
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Fig. 10. (Color online) Bell state pumping p i-> |^~)(T"|. Upper part: Pumping dynamics 
in Hilbert space, realized by two dissipative maps, under which two system qubits are 
irreversibly transferred from the +1 into the -1 eigenspaces of Z1Z2 and X1X2. Lower part: 
Schematics of the circuit decomposition of the two dissipative maps into unitary operations 
(i) - (iii), acting on the two system qubits S and an ancilla qubit playing the role of an 
environment E, followed by a dissipative reset (iv) of the ancilla. See main text for details. 
Figure adapted from Barreiro et al. (201 1) 

in the employed gate operations (see below). For pumping with unit probability 
(p = 1), the two qubits reach the target Bell state — regardless of their initial state 
— after only one cooling cycle, i.e., by a single application of each of the two 
maps. In contrast, in the limit p <^ I, the repeated application of this map generates 
dynamics according to a master equation (2) with Lindblad quantum jump operator 

C=iX2(l-ZiZ2). 



The map is implemented by a quantum circuit of three unitary operations (i)-(iii) 
and a dissipative step (iv), acting on two system qubits S and an ancilla which plays 
the role of the environment E (see lower part of Fig. 10): (i) Information about 
whether the system is in the +1 or -1 eigenspace of Z1Z2 is mapped by M(ZiZ2) 
onto the logical states |0) and |1) of the ancilla (initially in |1)): (ii) A controlled 
gate C(p) converts +1 into -1 eigenstates by flipping the state of the second qubit 
with probability p, where 

Cip) = mO\o®Ux,{p) + \l){l\o^l. 

Here, Uxiip) = exp(zQ'X2) and p = sin^ a controls the pumping probability, (iii) The 
initial mapping is inverted by A/"'(ZiZ2). At this stage, in general, the ancilla and 
system qubits are entangled, (iv) The ancilla is dissipatively reset to |1), which 
allows to carry away entropy to "cool" the two system qubits. The second map for 



24 



pumping into the -1 eigenspace of X1X2 is obtained from interchanging the roles of 
X and Z above. 

Experimental Bell state pumping - The described dynamics of "Bell state pumping" 
has been explored experimentally with three ions encoding the two system qubits 
and the ancilla qubit (Barreiro et al., 201 1) (see Fig. 2d). The unitary steps (i)-(iii) 
have been decomposed into the available set of coherent gate operations as shown 
in Fig. 2b and c. The dissipative reset of the ancilla qubit (iv) to its initial state |1) 
is realized by an addressed optical pumping technique, which leaves the quantum 
state of the system qubits unaffected (Schindler et al., 2011). The experimental 
results of various cycles of deterministic (p = 1) and probabilistic {p = 0.5) Bell 
state pumping are shown and discussed in Fig. 1 1 . 



2.3.2 Stabilizer Pumping and Ground State Cooling of the Toric Code Hamilto- 
nian 

The engineering of dissipative maps can be readily generalized to systems of more 
qubits. In particular, in the Rydberg simulator architecture (Weimer et al., 2010) a 
combination of coherent multi-atom Rydberg gates Ug (Eq. (4)) with optical pump- 
ing of ancillary control atoms allows one to implement collective dissipative many- 
particle dynamics in an a priori scalable system. As an example, we outline the 
engineering of dissipative dynamics for ground state cooling of Kitaev's toric code 
according to the plaquette and vertex four-body quantum jump operators given in 
Eq. (3). In direct analogy to the quantum circuit for Bell state pumping discussed in 
the previous section, four-qubit stabilizer pumping for a single plaquette is realized 
by a sequence of three unitary steps (shown in Fig. 12a), which are applied to the 
four system spins and the ancilla atom located at the center of the corresponding 
plaquette, followed by (iv) a dissipative reset of the ancilla qubit to its initial state. 

To this purpose, as for the simulation of coherent many -body dynamics (i) one first 
applies the mapping M (as specified in detail in the caption of Fig. 9) to coherently 
encode the information, whether the four system spins are in a -I- 1 or -1 eigenstate of 
the stabilizer Ap in the logical states of the auxiliary qubit, as schematically shown 
in Fig. 12b. (ii) Subsequently, a controlled spin flip onto one of the four system spins 
is applied, which converts a -1 ("high-energy") into a -i-l ("low-energy") eigenstate 
of Ap, with a certain, tunable probability determined by a phase (see Fig. 12a). 
(iii) After reversing the mapping M, the auxiliary qubit remains in the state 11)^., if 
one of the system spins has been flipped in the previous step (ii). Thus, (iv) finally 
addressed optical pumping resets the auxiliary ion from |l)c to its initial state |0)c, 
thereby guaranteeing that the auxiliary qubit factors out from the system dynamics 
and is "refreshed" for subsequent simulation steps. 

For small phases (f> (and thus small probabilities for pumping from the -1 into -l-l 
subspace of Ap in one step) and under a repeated application of this dissipative map. 
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Fig. 11. (Color online) Experimental Bell state pumping. Evolution of the Bell-state pop- 
ulations lO"^) (down triangles), |0~) (circles), |*F^) (squares) and |^^) (up triangles), a) 
Pumping process of an initially mixed state with probability p = I into the target Bell state 
\^~). Regardless of experimental imperfections, the target state population is preserved un- 
der the repeated application of further cooling cycles and reaches up to 91(1)% after 1.5 
cycles (ideally 100%). b) In a second experiment towards the simulation of master-equation 
dynamics, the probability is set at p = 0.5 to probe probabilistic cooling dynamics. In this 
case the target state is approached asymptotically. After cooling the system for 3 cycles with 
p = 0.5, up to 73(1)% of the initially mixed population cools into the target state (ideally 
88%). Error bars, not shown, are smaller than 2% (lo"). c) In order to completely character- 
ize the Bell-state cooling process, a quantum process tomography was performed (Nielsen 
and Chuang, 2000). As an example, the figure shows the reconstructed process matrix x 
(real part) for deterministic pumping with p = I, displayed in the Bell-state basis, describ- 
ing the deterministic pumping of the two ions after 1.5 cycles. The reconstructed process 
matrix has a Jamiolkowski process fidelity (Gilchrist et al., 2005) of 0.870(7) with the ideal 
dissipative process p i-> |*}'~){*P"|. The ideal process mapping any input state into the state 
|*F ) has as non-zero elements only the four transparent bars shown. Figure adapted from 
Barreiro et al. (2011). 

the density matrix p of the spin system evolves according to the Lindblad master 
Eq. (2) with the jump operators Cp given in Eq. (3) and the cooling rate k = (p^ It. 
Note, that the cooling also works for large phases 0; in fact, the most efficient dis- 
sipative state preparation is achieved with = nil, i.e. for deterministic pumping 
where an anyonic excitation, if it is present on the plaquette under consideration, 
is moved to a neighboring plaquette with unit probability. If this dynamics is ap- 
plied to all plaquettes of the lattice, it leads, as discussed above and illustrated in 
Fig. 7, to a dissipative random walk and pairwise annihilation of anyonic quasi- 
particle excitations, and thus in the long-time limit to a cooling of the system into 
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Fig. 12. (Color online) In a) Quantum circuit for the simulation of dissipative n-body inter- 
actions. In b) In a first step (i) the information whether the system spins reside in a + 1 or 
-I eigenstate of Ap is coherently mapped onto the logical states of the auxiliary qubit - in 
direct analogy to the simulation of coherent «-body interaction discussed above, (ii) Subse- 
quently a two-qubit gate Uz,ii4>) - |0)(Olc®l + U)(llc®^ with E = exp{i^cn) is applied. The 
"low-energy" +1 eigenstates of Ap are not affected by Uz,i as they have been mapped onto 
|0)c in step (i). In contrast - with probability p = sin^ (f> - the two-qubit gate induces a spin 
flip on the i-th system spin, if the system spins are in "high-energy" -1 eigenstates of Ap. 
(iii) The mapping M is reversed and (iv) finally, the auxiliary control qubit is incoherently 
reinitialized in state 10)^ by optical pumping. Controlling the angle in the quantum circuit 
allows one to realize either probabilistic cooling (0 «: k/2) described by a master equation 
with four-spin jump operators Cp as given in Eq. (3) or deterministic cooling (<p = n/2) as 
described by a discrete Kraus map of Eq. (6). 

the ground- state manifold. 

Many-body stabilizer pumping with trapped ions - Whereas for the described Ryd- 
berg simulator setup, all required components are not yet available in a single labo- 
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Fig. 13. (Color online) Experimental four-qubit stabilizer pumping, which can be re- 
garded as dissipative ground state preparation of one plaquette of Kitaev's toric code 
(Kitaev, 2003). a) Schematic of the four system qubits to be cooled into the GHZ state 
(10000) + 11111))/ yl, which is uniquely characterized as the simultaneous eigenstate with 
eigenvalue +1 of the shown stabilizers, b) Reconstructed density matrices (real part) of 
the initial mixed state pmixed and subsequent states pi, 2,3,4 after sequentially pumping the 
stabilizers Z1Z2, Z2Z3, Z3Z4 and X1X2X3X4. Populations in the initial mixed state with 
qubits / and j antiparallel, or in the -1 eigenspace of the ZiZj stabilizer, disappear after 
pumping this stabilizer into the + 1 eigenspace. For example, populations in dark blue dis- 
appear after ZiZ2-stabilizer pumping. A final pumping of the stabilizer X1X2X3X4 builds 
up the coherence between |0000) and |1111), shown as red bars in the density matrix 
of p4. The reconstructed density matrices for the initial and subsequent states arising 
in each step have a fidelity, or state overlap (Jozsa, 1994), with the expected states of 
{79(2),89(1),79.7(7),70.0(7),55.8(4))%. c) Measured expectation values of the relevant sta- 
bilizers; ideally, non-zero expectation values have a value of -i-l. d) Evolution of the mea- 
sured expectation values of the relevant stabilizers for repetitively pumping an initial state 
|1111) with probability p = 0.5 into the -1 eigenspace of the stabilizer X1X2X3X4. The in- 
cremental cooling is evident by the red line fitted to the pumped stabilizer expectation value. 
The evolution of the expectation value (X1X2X3X4) for deterministic cooling (p = 1) is also 
shown. The observed decay of (ZiZj) is due to imperfections and detrimental to the pump- 
ing process. Error bars in (c) and (d), ±lcr. Figure reprinted with permission from Barreiro 
et al. (2011). Copyright 2011 by MacMillan. 

ratory, Barreiro et al. (201 1) demonstrated the described four-qubit stabilizer pump- 
ing in a proof-of-principle experiment with 5 trapped ions. Specifically, pumping 
dynamics into a four-qubit Greenberger-Horne-Zeilinger (GHZ) state (|0000) -l- 
|1 1 1 1))/ V2 was realized. This state can be regarded as the ground state of a min- 
imal instance of Kitaev's toric code, consisting of a single square plaquette, as 
sketched in Fig. 13a. The state is uniquely characterized as the simultaneous eigen- 
state of the four stabilizers Z1Z2, Z2Z3, Z3Z4 and X1X2X3X4, all with eigenvalue -l-l. 
Therefore, cooling dynamics into the GHZ state are realized by four consecutive 
dissipative steps, each pumping the system into the -l-l eigenspaces of the four sta- 
bilizers (Fig. 13b-d). In a system of 4-1-1 ions encoding the four system spins and an 
ancillary qubit, such cooling dynamics has been realized in analogy with the Bell- 
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state pumping discussed in Sect. 2.3. 1 . Here, however, the experimental complexity 
is considerably larger, as the circuit decomposition of one cooling cycle involves 
16 five-ion entangling M0lmer-S0rensen gates, 20 collective and 34 single-qubit 
rotations; further details in Barreiro et al. (201 1). 

Open- vs. closed-loop control scenarios - In the discussed examples of engineering 
of dissipative dynamics for Bell-state and four-qubit stabilizer pumping the avail- 
able quantum resources were used by coupling the system qubits to an ancilla qubit 
by a universal set of gates. Such set was constituted by entangling multi-ion MS 
gates in combination with single-ion rotations (Barreiro et al., 2011), or the the 
mesoscopic Rydberg gate (Miiller et al., 2009) in combination with single-atom 
gates in the Rydberg simulator architecture (Weimer et al., 2010). The engineered 
environment was here represented by ancilla ions or Rydberg atoms, undergoing 
optical pumping by dissipative coupling to the vacuum modes of the radiation 
field. Note that in the described scenario, the ancilla qubit remains unobserved, 
representing an open-loop dynamics. For such open quantum systems, though, it 
was noted in Bacon et al. (2001); Lloyd and Viola (2001) that using a single an- 
cilla qubit the most general Markovian open-system dynamics cannot be obtained 
with a finite set of non-unitary open-loop transformations. However, such a uni- 
versal dynamical control can be achieved through repeated application of coherent 
control operations and measurement of the auxiliary qubit, followed by classical 
feedback operations onto the system. In the trapped-ion experiments in Barreiro 
et al. (2011) the simulation toolbox was complemented by the demonstration of 
a quantum-non-demolition (QND) measurement of a four-qubit stabilizer operator 
via an auxiliary qubit. In combination with classical feedback (Riebe et al., 2008), 
such QND readout operations provide the basis for such closed-loop dynamics. 

Furthermore, in the context of quantum error correction, QND measurements of 
stabilizer operators constitute a crucial ingredient for the realization of quantum 
error-correcting codes (Steane, 1996; Calderbank and Shor, 1996). Such readout 
operations correspond to error syndrome measurements, and the obtained infor- 
mation can be classically processed and used to detect and correct errors (Dennis 
et al., 2002). For instance, in Miiller et al. (2011) it is explicitly worked out how 
minimal instances of complete topogical quantum error correcting codes (Bombin 
and Martin-Delgado, 2006) can be realized with the currently available toolbox for 
open-system quantum simulation with trapped ions (Barreiro et al., 201 1). 



2.3.3 Digital Simulation of a U(l) Lattice Gauge Theory 

The above analysis for the coherent simulation and ground state cooling of Kitaev's 
toric code can be extended to a large class of interesting models. In Weimer et al. 
(2011) it is discussed how the digital Rydberg simulator architecture enables the 
simulation of Heisenberg-like spin models, and in principle also fermionic Hubbard 
models, by mapping lattice fermions to a spin Hamiltonian involving many-body 
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interactions that can be realized in the Rydberg simulator. 

Three-dimensional U(l) lattice gauge theory - The toric code is the ground state 
of the frustration-free, exactly solvable toric code Hamiltonian involving four-qubit 
plaquette and vertex interactions (Kitaev, 2003). It belongs to the class of stabilizer 
states and exhibits Abelian topological order. It is also possible to provide (digi- 
tal) simulation protocols for the simulation of coherent many-body dynamics and 
ground state preparation of more complex spin models. In Weimer et al. (2010) 
such a protocol was developed for the example of a three-dimensional [/(l)-lattice 
gauge theory (Kogut, 1979) and it was shown how to achieve dissipative ground 
state preparation also for such a complex system. Such models have attracted inter- 
est in the search for 'exotic' phases and spin liquids (Moessner and Sondhi, 2001; 
Motrunich and Senthil, 2002; Hermele et al., 2004; Levin and Wen, 2005a,b). 

Fig. 14a shows the setup of the U(l) lattice gauge theory. Spins are located on the 
edges of a three-dimensional cubic lattice and interact via the many-body Hamilto- 
nian 

H = uJ]{Slf -JJ]b„ + VN^k. (8) 

o p 

The first term with Sg = ^jteo (^[ ^ describes pairwise two-body interactions of six 
spins located at the corners of octahedra, located around the vertices of the square 
lattice (see the spins connected by red lines in Fig. 14a). The inequality U » \V\, \J\ 
defines a low-energy sector of the theory, which consists of spin configurations with 
an equal number of three up and three down spins, i.e., states with vanishing total 
spin 5^ on each octahedron. The second term describes a ring-exchange interac- 
tion Bp = S'^S2S'^S^ -I- 5 [5 2 5 3^4 of four spins located around each plaquette of 
the lattice (see green plaquette in Fig. 14a); here Sf = (cr^ + icr])/2. This inter- 
action flips the state of four plaquette spins with alternating spin orientation, e.g., 
Bp\OlOl)p = |1010)p, and leaves other states unchanged, e.g., 5p|1001)p = 0. Note 
that while the ring-exchange interaction term commutes with the S I spin constraint 
terms, ring-exchange terms on neighboring plaquettes do not commute. 

The last term of the Hamiltonian of Eq. (8) counts the total number of flippable 
plaquettes A^rk = T^pB],. It is introduced since at the so-called Rokhsar-Kivelson 
point with J = V, the system becomes exactly solvable (Rokhsar and Kivelson, 
1988). If one identifies each spin up with a "dimer" on a link of the lattice, all 
states satisfying the low-energy constraint of vanishing 5 ^ on all octahedra can be 
viewed as an "allowed" dimer covering with three dimers meeting at each site of 
the cubic lattice. Fig. 14b shows how the Bp ring exchange interaction term flips 
one dimer covering into another. Within this dimer description, the ground state 
at the Rokhsar-Kivelson point is given by the condensation of the dimer coverings 
(Levin and Wen, 2005b), i.e., the equal weight superposition of all allowed dimer 
coverings. It has been suggested that in the non-solvable parameter regime < V < 
J of interest the ground state of the system is determined by a spin liquid smoothly 
connected to the Rokhsar-Kivelson point (Hermele et al., 2004). 
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a Bp = S+S2 S+S^ +h.c 




Fig. 14. (Color online) Simulation of a three-dimensional U(l) lattice gauge theory, a) 
Spins located on the edges of a cubic lattice interact via a six-spin low-energy constraint 
term S^ (indicated by red links), which imposes the condition of an equal number of three 
up- and three down-spins on each octahedron, and via a four-spin ring-exchange plaquette 
interaction Bp (green links) - see Eq. (8) and details in the main text, b) In the language of 
dimer coverings, the ring exchange terms Bp coherently convert flippable plaquettes from 
one configuration into another, c) Numerical simulation of the cooling into the ground state 
at the Rokhsar-Kivelson point V = 7 for a system of 4 unit cells (12 spins): The cooling into 
the low-energy subspace defined by the octahedra constraints can be realized in analogy to 
the cooling of the toric code (Weimer et al., 2010); alternatively one can directly start 
the protocol in an initial (classical) state, which satisfies all constraints. The inset shows 
the cooling into the equal-weight superposition of all dimer coverings starting from an 
initial state which already satisfies the S ^-constraints, d) Coherent time evolution from the 
Rokhsar-Kivelson point with a linear ramp of the Rokhsar-Kivelson term V{t) = (l-tJ/lO): 
the solid line denotes the exact ground state energy, while dots represent the digital time 
evolution during an adiabatic ramp for different phases during each time step (Weimer 
et al., 2010). The differences arise from Trotter expansion errors due to non-commutative 
terms in the Hamiltonian (8). Parts (b) - (d) of the figure adapted from Weimer et al. (2010)] 



Simulation protocol - To reach the < V < 7 phase of interest, the idea is to 
(i) implement dissipative dynamics, which first cools the system at the Rokshar- 
Kivelson point (7 = V) into the ground state given by the symmetric superposition 
of dimer coverings, and (ii) subsequently to slowly decrease the strength of the 
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Rokshar-Kivelson term V Nrk in the Hamiltonian (8) such that the ground state is 
adiabatically transformed into the quantum phase of interest: 

(i) If one starts in some initial state, which satisfies the SI on all octahedra, the 
condensation of the dimer coverings can be achieved by dissipative dynamics ac- 
cording to plaquette jump operators 

Cp=^-CT'^[l-B^]Bp. (9) 

The jump operator Cp has by construction two dark states, which are the and 
+ 1 eigenstates of Bp. The eigenstates correspond to a non-flippable plaquette 
(e.g. Cp|1001)p = 0), while the +1 eigenstate is the equal- weight superposition of 
the original dimer covering and the dimer covering obtained by flipping the plaque- 
tte, |1010)p -I- |0101)p. Finally, the jump operator Cp transforms the third eigenstate 
with eigenvalue -1 into the -i-l eigenstate. As a consequence, as Fig. 14c illustrates, 
under this dynamics acting on all plaquettes of the cubic lattice, for long times the 
system asymptotically approaches the ground state consisting of the symmetric su- 
perposition of all allowed dimer coverings. 

(ii) Subsequently, this ground state is transformed adiabatically into the phase at < 
y < 7 by slowly ramping down the Rokshar-Kivelson term. Such adiabatic passage 
can be realized by decomposing the coherent dynamics according to the Hamilto- 
nian with the time-dependent Rokshar-Kivelson term V(t) Nr^ into small Trotter 
time steps (conceptually similar to the simulation of two-spin time-dependent Trot- 
ter dynamics discussed in Sect. 2.1.2). The different curves in Fig. 14d indicate 
deviations of the simulated adiabatic passage from the exact dynamics due to Trot- 
ter errors originating from the non-commutativity of terms in the Hamiltonian (8). 

The Hamiltonian terms (8) and quantum jump operators (9) for the simulation of 
the U(l) lattice gauge theory are more complex than the ones for ground state cool- 
ing and Hamiltonian dynamics according to the toric code Hamiltonian. However, 
in the Rydberg simulator architecture they can also be implemented by combina- 
tions of many-atom Rydberg gates and optical pumping of ancilla qubits, which are 
located on the plaquettes and corners of the qubit lattice; see Weimer et al. (2010) 
for details and explicit circuit decompositions. 



2.4 The Effect of Gate Imperfections on Digital Quantum Simulation 



Imperfect gate operations in the quantum circuits which are used to implement co- 
herent and dissipative steps of time evolution according to discrete Kraus maps (6) 
or many-body master equations (1) lead to deviations of the actually realized sys- 
tem from the envisioned dynamics. In the simulation of many-body dynamics for a 
given time t via a Trotter decomposition this leads in practice to a trade-off": On the 
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one hand, the number of simulation steps n according to small time intervals tin 
should be chosen large, in order to keep the effect of Trotterization errors originat- 
ing from non-commuting terms small. On the other hand, the practical implemen- 
tation of each time step has a certain cost in terms of resources and is associated to 
a certain experimental error, which favors the decomposition of the simulated time 
dynamics into a not too large number of steps. 

Small imperfections typically provide in leading order small perturbations for the 
simulated Hamiltonian dynamics and weak additional dissipative terms. The spe- 
cific form is strongly dependent on the particular implementation platform and its 
dominant error sources; see the analysis in Diir et al. (2008) for a general discus- 
sion. For the Rydberg quantum simulator architecture (Weimer et al., 2010) the 
influence of errors in the multi-atom Rydberg gate (Miiller et al., 2009) on the sim- 
ulation of Kitaev's toric code Hamiltonian and ground state cooling in this model 
has been analyzed: Fig. 7b shows that in the presence of small gate imperfections 
the desired cooling into the ground state of the model is accompanied by weak, 
unwanted heating processes, such that in the long-time limit a finite anyon density 
remains present in the many-body system. Such eff"ects have also been observed 
experimentally in the dissipative state preparation of a minimal system of one pla- 
quette of the toric code with trapped ions (Barreiro et al., 2011), as discussed in 
Sect. 2.3.2: Fig. 13d shows that under repeated pumping into the -1 eigenspace of 
the four-qubit- stabilizer X1X2X3X4, the expectation value of the two-qubit stabiliz- 
ers ZiZj, which should ideally be unaffected by the XiX2X3X4-pumping and should 
remain at their initial value of -i-l, undergo a decay. This detrimental effect can be 
interpreted as "heating processes" due to experimental imperfections in the under- 
lying quantum circuits; see also Miiller et al. (2011) where a theoretical model- 
ing of these errors is discussed. However, the thermodynamic properties (quantum 
phases) and dynamical behaviour of a strongly interacting many-body system are 
in general robust to small perturbations in the Hamiltonian; e.g., the stability of the 
toric code for small magnetic fields has recently been demonstrated (Vidal et al., 
2009). Consequently, small imperfections in the implementation of the gate oper- 
ations leading to deviations from the ideal simulated dynamics are expected to be 
tolerable. 



3 Engineered Open Systems with Cold Atoms 



As anticipated in Sect. 1, here we will be interested in a scenario where many-body 
ensembles of cold atoms are properly viewed as open quantum systems, in a setting 
familiar from quantum optics: A system of interest is coupled to an environment, 
giving rise to dissipative processes, and is additionally driven by external coherent 
fields. This creates a non-equilibrium many-body setting without immediate coun- 
terpart in condensed matter systems. In particular, in the first part of this section, 
we point out how the conspiracy of laser drive and dissipation can give rise to off- 
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diagonal long-range order, a trademark of macroscopic quantum phenomena. We 
also argue how this can be achieved via proper reservoir engineering, in this way 
fully extending the notion of quantum state engineering from the Hamiltonian to 
the more general Liouvillian setting, where controlled dissipation is included. 

In the following parts of this section, we will give accounts for further central as- 
pects of this general setting. In the context of atomic bosons, we point out in which 
sense these systems indeed constitute a novel class of artificial out-of-equilibrium 
many-body systems, by analyzing a stationary state phase diagram resulting from 
competing unitary and dissipative dynamics. In the context of atomic fermions, 
we present a dissipative pairing mechanism which builds on a conspiracy of Pauli 
blocking and dissipative phase locking, based on which we argue that such sys- 
tems may provide an attractive route towards quantum simulation of important 
condensed matter models, such as the Fermi-Hubbard model. We then explore the 
possibilities of dissipatively realizing topological phases in the lab, and elaborate 
on the specific many-body properties of such dissipatively stabilized states of mat- 
ter. 

The results presented here highlight the fact that the stationary states of such 
driven-dissipative ensembles, representing flux equilibrium states far from ther- 
modynamic equilibrium, feature interesting many -body aspects. This places these 
systems in strong contrast to the dynamical non-equilibrium phenomena which are 
currently actively investigated in closed systems in the cold atom context, focusing 
on thermalization (Gasenzer et al., 2005; Cramer et al., 2008; Rigol et al., 2008; 
Kinoshita et al., 2006; Hofi'erberth et al., 2007; Trotzky et al., 2011) and quench 
dynamics (Calabrese and Cardy, 2006; KoUath et al., 2007; Greiner et al., 2002b; 
Sadler et al., 2006). 



3.1 Long-Range Order via Dissipation 



3.1.1 Driven-Dissipative BEC 

Qualitative picture: Dark states in single- and many-particle systems - For long 
times, a system density matrix governed by Eq. (1) will approach a flux equilib- 
rium stationary state, p{t) -^ pss, in the presence of dissipation, which generically 
is a mixed state. However, under suitable circumstances the stationary state can 
be a pure state, p^s = \D) {D\. In the language of quantum optics, such states \D) 
are called dark states. A familiar example on the level of single particles is optical 
pumping or dark state laser cooling to subrecoil temperatures (Aspect et al., 1988; 
Kasevich and Chu, 1992), illustrated in Fig. 15a: By coherently coupling two de- 
generate levels to an auxiliary excited state with antisymmetric Rabi frequencies 
+Q., from which spontaneous emission leading back to the ground states occurs 
symmetrically, a dark state is given by the symmetric superposition of the ground 
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states. For sufficient detuning, it is then clear that the population will entirely end 
up in this dark state decoupled from the light field. In our setting, we replace the 
internal degrees of freedom of an atom by external, motional degrees of freedom, 
realized schematically by an optical potential configuration with an intermediate 
site on the link between degenerate ground states, cf. Fig. 15b. Below we will dis- 
cuss how to realize the relevant driving and decay processes. Clearly, the same ar- 
guments then lead to a phase locked, symmetric superposition dark state as above, 
i.e. {a\ -\- fl2)|vac) in a second quantization notation. However, two generalizations 
follow immediately: First, the levels (lattice sites) can be populated with bosonic 
degrees of freedom, i.e. there is no limit on the occupation number. Second, and 
most natural in an optical lattice context, the "dark state unit cell" can be cloned 
in a translation invariant way to give a complete lattice setting, in one or higher 
dimensions. The key ingredient is antisymmetric drive of each pair of sites, and the 
spontaneous decay back to the lower states, as depicted in Fig. 15c. The phase is 
then locked on each two adjacent sites, such that eventually only the symmetric su- 
perposition over the whole lattice persists. This is the only state not being recycled 
into the dissipative evolution. This state is nothing but a Bose-Einstein condensate 
(BEC) with a fixed but arbitrary particle number A'^, which for a one-dimensional 
geometry with M sites depicted in Fig. 15c reads 



i 

In consequence, quantum mechanical long range order is built up from quasilocal, 
number conserving dissipative operations. The system density matrix is purified, in 
that a zero entropy state is reached from an arbitrary initial density matrix, as will 
be discussed next. 

Driven-dissipative BEC as unique stationary state - Here we make the above intu- 
itive picture more precise by discussing the Lindblad jump operators which drive 
into the BEC state Eq. (10). In a slight generalization, we consider the dynamics 
of A'^ bosonic atoms on a J-dimensional lattice with spacing a and M'^ lattice sites, 
and lattice vectors e^. For simplicity, we first address the purely dissipative case of 
Eq. (1) and set H = 0. The goal is then achieved by choosing the jump operators 
Eq. (2) as 

C/J = <^U = {^1 + '^y) {^i - '^j) ' (11) 

acting between each pair of adjacent lattice sites /3 = (i, j) with an overall dissipa- 
tive rate k^ = Kjj = k. Because the annihilation part of the normal ordered operators 
Cjs commute with the generator of the BEC state 2; ^J^ we have 



(fl,-a,)|BEC) = OV</,j), (12) 

making this state indeed a many-body dark state (or dissipative zero mode) of the 
Liouville operator defined with jump operators Eq. (11). 
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\D) ~ \g+i) + \g-i) \B) ~ \g+i) - \g_i) 
dark state bright state 



b) ^\j^ 




{a\ + al) |vac) (a| - al) |vac) 
symmetric anti-symmetric 







(;/^ Hilbert space 




dark subspace 



Fig. 15. (Color online) Dark states in many-body systems from an analogy with opti- 
cal pumping: a) A coherently driven and spontaneously decaying atomic A-system with 
metastable excited state has the symmetric superposition of the degenerate ground states as 
a dark state for antisymmetric driving, b) The internal degrees of freedom are replaced by 
external degrees of freedom, such as the sites of an optical superlattice, with the same con- 
sequences once antisymmetric driving and spontaneous emission are properly engineered, 
c) The unit cell is naturally cloned in a translation invariant lattice setting. The symmetric 
phase locking on each pair of sites generates coherence over the whole system, correspond- 
ing to a fixed number BEC. d) Uniqueness: If the dark subspace consists of one dark state 
only, and no subspace exists which is left invariant under the set of jump operators, the 
many -body density matrix converges to the dark state irrespective of its initial condition. 

From the explicit form of the jump operators, we see that the key for obtaining a 
dark state with long range order is a coupling to the bath which involves a current or 
discrete gradient operator between two adjacent lattice sites. The temporally local 
jump operator Cjj describes a pumping process, where the annihilation part a, - cij 
removes any anti-symmetric (out-phase) superposition on each pair of sites (/, j), 
while aj -I- a^. recycles the atoms into the symmetric (in-phase) state. As anticipated 
above, this process can thus be interpreted as a dissipative locking of the atomic 
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phases of every two adjacent lattice sites, in turn resulting into a global phase lock- 
ing characteristic of a condensate. 

We also note from Eq. (12) that the dark state property of |BEC) is mainly deter- 
mined by the annihilation part of the jump operator. In fact, any linear combination 
of a],a^- of recycling operators will work, except for a hermitian Cij, i.e. for the 

combination a] -a}-. In this case, the dissipative dynamics would result in dephas- 
ing instead of pumping into the dark state. This case is then qualitatively similar 
to the generic situation in atomic physics. There, a bath typically couples to the 
atomic density with jump operators n, = aja,, as in the case of decoherence due to 
spontaneous emission in an optical lattice (Pichler et al., 2010), or for coUisional 
interactions. 

We now discuss the uniqueness of the stationary dark state. The following two re- 
quirements have to be fulfilled to ensure uniqueness (in the absence of Hamiltonian 
dynamics) (Baumgartner et al., 2007; Kraus et al., 2008): 

(i) The dark subspace is one-dimensional, i.e. there is exactly one normalized dark 
state |D), for which 

C^|D)=OV;0. (13) 

(ii) No stationary solutions other than the dark state exist. 

In the above example, so far we have only argued that the EEC state is a dark 
state. However, it is easily seen that no other dark states are present, since the non- 
hermitean creation and annihilation operators can only have eigenvalue zero on an 
//-particle Hilbert space. In particular, the creation part a'^ -I- a}- never has a zero 
eigenvalue, as it acts on an A^ - 1 particle Hilbert space. We can therefore focus 
on the annihilation part alone, where the Fourier transform 2^(1 ~ e"''^')^q reveals 
indeed exactly one zero mode at q = 0. As to (ii), uniqueness of the dark state as 
a stationary state is guaranteed if there is no other subspace of the system Hilbert 
space which is left invariant under the action of the operators c^ (Baumgartner 
et al., 2007; Kraus et al., 2008). This can be shown explicitly for the example above 
(Kraus et al., 2008). More generally, it can be proved that for any given pure state 
there will be a master equation so that this state becomes the unique stationary 
state. Uniqueness is a key property: under this circumstance, the system will be 
attracted to the dark state for arbitrary initial density matrix, as illustrated in Fig. 
15d. These statements remain true for a Hamiltonian dynamics that is compatible 
with the Lindblad dynamics, in the sense of the dark state being an eigenstate of 
the Hamiltonian, }i\D) = E \D). One example is the addition of a purely kinetic 
Hamiltonian, since //o |BEC) = A/^6q=o |BEC), where 6q = 272^sin^qei/2 is the 
single particle Bloch energy for quasimomentum q. 

Finally, we remark that as a consequence of the symmetry of global phase rotations 
exerted by e'^^ on the set of jump operators (i.e. [cp, A'^] = V/3, where N = 2,- a^a,- 
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superfluid 
' reservoir. 



Fig. 16. (Color online) Cold atom implementation of a driven-dissipative condensate: a) 
A coherently driven lattice gas is immersed in a surrounding condensate, b) Schematic 
realization of the effective dissipative process in an optical superlattice, which provides for 
excited states gapped by s and localized on the links of neighboring lattice sites {ij): A 
Raman laser couples the ground- and excited bands with effective Rabi-frequency O and 
detuning ^ - u - e from the inter-band transition. Only the antisymmetric component of 
atoms on neighboring lattice sites is excited to the upper band due to the spatial modulation 
of the Raman-laser. The inter-band decay with a rate Y back to the lower band is obtained 
via the emission of Bogoliubov quasiparticle excitations into the surrounding BEC. Figure 
adapted from Diehl et al. (2008). 

is the total particle number operator), which is present microscopically, any break- 
ing of this symmetry must occur spontaneously. This gives room for concept of 
spontaneous symmetry breaking to be applicable in the thermodynamic limit for 
such driven-dissipative systems. 



3.1.2 Implementation with Cold Atoms 

Before sketching an explicit implementation scheme of the above dynamics, we 
point out that the existence of a microscopic scale, where a description of the sys- 
tem in terms of a temporally local evolution equation is possible, is far from obvious 
in a many-body context. In fact, in usual condensed matter settings, typical baths 
have arbitrarily low energies which can be exchanged with a given many-body 
system of interest, giving rise to temporally non-local memory kernels in the de- 
scription of environmental effects. Instead, the validity of the master equation rests 
on the Bom-Markov approximation with system-bath coupling in rotating wave ap- 
proximation. This means that the bath is gapped in a condensed matter language. 
For typical quantum optics settings, these approximations are excellent because the 
(optical) system frequencies providing for the gap are much larger than the decay 
rates. Below we argue how to mimic such a situation in an optical lattice context. At 
the same time, this setting makes clear the need for external driving in order to pro- 
vide the energy necessary to access the decaying energy levels. The validity of this 
combination of approximations then fully extends the scope of microscopic control 
in cold atom systems from unitary to combined unitary-dissipative dynamics. 
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A concrete possible implementation in systems of cold bosonic atoms a builds 
on the immersion of a coherently driven optical lattice system into a large BEC 
of atoms b (Griessner et al., 2006), cf. Fig. 16a. In order to realize the key A- 
configuration, we consider a superlattice setting as illustrated in one-dimensional 
geometry in Fig. 15a, with an additional auxiliary lattice site on each of the links. 
The optical lattice corresponding to a single link is shown in Fig. 16b, where the 
A-system is implemented with the two Wannier functions of lattice sites 1 and 2 
representing two ground states, and the auxiliary state in the middle representing 
an excited state. In order to achieve the annihilation part of the jump operator, we 
drive this three-level system by Raman transitions from the two ground to the ex- 
cited states with Rabi frequencies Q. and -Q, respectively. This could be realized 
in a translation invariant way for the whole lattice by, e.g., a commensurate ratio 
of lattice and Raman laser wavelengths. Ataman = 2/liatt, which would guarantee the 
relative sign via a ;r-phase shift for the Rabi frequency. In the next step, the dissipa- 
tion needs to be introduced. To this end, the coherently driven system is placed into 
a large BEC reservoir. This condensate interacts in the form of a conventional s- 
wave contact potential with interspecies scattering length aab with the lattice atoms 
a, and acts as a bath of Bogoliubov excitations. Such a coupling provides an effi- 
cient mechanism for decay of atoms a from the excited to the lower Bloch band 
by emission of Bogoliubov quasiparticles. This replaces photon emission in a con- 
ventional quantum optics situation. The conspiracy of coherent drive and dissipa- 
tion explained here also gives rise to the physical picture of the coherence of the 
driving laser beam being imprinted onto the matter system - any deviation from 
the above commensurability condition would be reflected in a length scale in the 
driven-dissipative BEC. We note however, that the ratio of wavelengths can be con- 
trolled with high precision in state-of-the-art experiments. 

In the presence of a large condensate, linearization of the system-bath interaction 
around the bath condensate expectation value, together with the harmonic bath of 
Bogoliubov excitations, realizes the generic system-bath setting of quantum optics. 
In particular, a key element is the presence of the largest energy scale provided 
by the Hubbard band separation e (cf. Fig. 16b), ensuring the validity of Born- 
Markov and rotating wave approximations. This in turn leads to a temporally local 
master equation description. As long as this scale exceeds the bath temperature 
6 » Tbec, the occupation of modes at these energies is negligible and the BEC 
thus acts as an effective zero temperature reservoir. At the same time, the role of 
coherent driving with energy oj in order to bridge the energy separation of the two 
bands becomes apparent. The fact that energy is constantly pumped into the system 
in our driven-dissipative non-equilibrium setting highlights the fact that our setting 
can indeed realize states of zero entropy, or in practice an entropy substantially 
lower than the surrounding reservoir gas, without conflicting with the second law 
of thermodynamics. 

If we further specialize to the limit of weak driving Q «; A, where A = (x> - 6 is a 
detuning from the upper Hubbard band, adiabatic elimination of the excited Bloch 
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band results in a master equation generated by jump operators of the type (11). 
In this case, on the full lattice, the laser excitation to the upper band ~ (a, - Gj) 
for each pair of sites is followed by immediate return of the atoms into the lowest 
band, which generically happens in a symmetric fashion such as ~ {a] + aX in 
this way realizing jump operators of the form of Eq. (11). Details of the return 
process, however, depend on the Bogoliubov excitation wavelength in the bath: 
For wavelength A\, larger or smaller than the optical lattice spacing a, spontaneous 
emission is either correlated or uncorrected. However complicated, the existence 
of a dark state in the present case is guaranteed by (a, - aj) |BEC) = 0, a property 
which follows from the laser excitation step alone. We will therefore concentrate 
below on the jump operators defined in Eq. (11). 

Finally, we emphasize that the basic concept for the dissipative generation of long- 
range order in many-particle systems can be explored in very different physical 
platforms beyond the cold atom context, offering additional opportunities for im- 
plementations. For example, microcavity arrays have been identified as promising 
candidates for the realization of the above dynamics with state-of-the-art technol- 
ogy (Marcos et al., 2012), where the bosonic degrees of freedom are realized by 
microwave cavity photons. The auxiliary system is there realized by two interacting 
superconducting qubits, which are placed between two neighboring microwave res- 
onators. The symmetric and antisymmetric superposition modes of the resonators 
are coupled to the qubit system and the dissipative step is realized naturally via 
spontaneous decay of the latter. 

In an even broader context, also different kinds of intrinsically quantum mechan- 
ical correlations, such as entanglement, can be targeted dissipatively. Examples 
have been discussed in trapped ion systems above. In addition, in a recent break- 
through experiment entanglement has been generated dissipatively between two 
macroscopic spin ensembles (Krauter et al., 2011; Muschik et al., 2011). On the 
theory side, creation of atomic entanglement has been proposed in the context of 
optical cavities (Kastoryano et al., 201 1), and the generation of squeezed states of 
matter has been investigated for the case of macroscopic two-mode boson ensem- 
bles (Makela and Watanabe, 201 1). Furthermore, dissipation has been proposed as 
a means to purify many-body Fock states as defect- free registers for quantum com- 
puting with cold atoms (Pupillo et al., 2004; Brennen et al., 2005), as well as to en- 
force three-body constraints in Hamiltonian dynamics (Daley et al., 2009; Kantian 
et al., 2009; Diehl et al., 2010a; Roncaglia et al., 2010). Recent landmark experi- 
ments have used it to build strong correlations in, and thus to stabilize, a metastable 
weakly interacting molecular gas in one-dimension (Syassen et al., 2008; Porto, 
2008). 

So far we have discussed the proof- of-principle for the concept of state engineering 
in many-particle systems by tailored dissipation in the conceptually simplest exam- 
ple, the driven-dissipative BEC. In the following subsections, we will review differ- 
ent research directions which address many-body aspects in such systems, where 
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dissipation acts as a dominant resource of dynamics: In the context of bosonic sys- 
tems, we present a dynamical phase transition resulting from the competition of 
the engineered Liouville- with a Hamiltonian dynamics, defining a novel class of 
interacting non-equilibrium many-body systems with interesting stationary states. 
The phase transition is seen to share features of both quantum and classical phase 
transitions, and we identify an intriguing phase where global phase rotation and 
translation symmetry are simultaneously broken spontaneously. In the context of 
atomic fermions, we discuss a dissipative pairing mechanism, which is operative 
in the absence of attractive forces and allows us to target states of arbitrary sym- 
metry, such as d-wave paired states in two dimensions. Beyond the identification 
of this new far-from-equilibrium pairing mechanism, this makes dissipative state 
engineering potentially relevant for the experimental efi'orts towards the quantum 
simulation of the two-dimensional Fermi-Hubbard model, where the ground state 
is believed to have pairing with d-wave symmetry away from half filling. Finally, 
we show how dissipation engineering can be used in order to reach fermionic states 
with topological order dissipatively. While so far topological phases have been ex- 
clusively discussed in a Hamiltonian context, we develop here a dissipative coun- 
terpart for such phases. We discuss the associated phenomena resulting when such 
systems are suitably constrained in space, such as the emergence of unpaired Ma- 
jorana edge modes. 



3.2 Competition of Unitary and Dissipative Dynamics in Bosonic Systems 



Motivation - In a Hamiltonian ground state context, a quantum phase transition 
results from the competition of two non-commuting parts of a microscopic Hamil- 
tonian H = Hi + gH2, if the ground states for g <^ I and g » 1 have diff"erent 
symmetries (Sachdev, 1999). A critical value gc then separates two distinct quan- 
tum phases described by pure states, while in thermodynamic equilibrium for finite 
temperature this defines a quantum critical region around gcinaT vs. g phase dia- 
gram. Classical phase transitions may occur for fixed parameter g by increasing the 
temperature, and can be viewed as resulting from the competition of the specific 
ground state stabilized by the Hamiltonian vs. the completely mixed structureless 
infinite temperature state. In contrast, here we study a non-equilibrium situation, in 
which there is a competition between a Hamiltonian and a dissipative dynamics. 
We extract the complete steady state phase diagram, revealing that the resulting 
transitions share features of quantum phase transitions, in that they are interac- 
tion driven, and classical ones, in that the ordered phase directly terminates into 
a strongly mixed state. It contains an extended region where global phase rotation 
and translation symmetry are both broken spontaneously, as a consequence of a 
subtle renormalization eff"ect on the complex excitation spectrum of the low-lying 
modes. In addition, we study the dynamical critical behavior in the long-time limit 
of the combined unitary and dissipative evolution. 
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Those aspects underpin the fact that the driven-dissipative systems investigated 
here add a new class of non-equilibrium stationary states to those which have 
been studied so far. One prominent example is certainly electron systems in con- 
densed matter, which are exposed to a bias voltage (Kamenev and Levchenko, 
2009). In this context, also characteristic many-body behavior such as the effect of 
non-equilibrium conditions on quantum critical points has been investigated (Mitra 
et al., 2006). Further routes of driving many-body systems out of thermodynamic 
equilibrium are discussed in the context of exciton-polariton Bose-Einstein con- 
densates (Moskalenko and Snoke, 2000; Kasprzak et al., 2006), or more recently 
in driven noisy systems of trapped ions or dipolar atomic gases (Dalla Torre et al., 
2010,2011). 



3.2.1 Dynamical Phase Transition 

Model and Analogy to Equilibrium Quantum Phase Transition - We now extend 
the purely dissipative dynamics leading to a BEC state determined by Eq. (1 1) by 
the generic Hamiltonian in optical lattice systems, the Bose-Hubbard Hamiltonian: 

d,p = -i[H,p]+£[p], (14) 

This Hamiltonian is defined with the parameters 7, the hopping amplitude, and [/, 
the onsite interaction strength; h( = b\bi is the number operator for site £. Its ground 
state physics provides a seminal example for a quantum phase transition in the cold 
atom context (Fisher et al., 1989; Jaksch et al., 1998b; Greiner et al., 2002a; Bloch 
et al., 2008): For a given chemical potential yu, which in equilibrium fixes the mean 
particle density n, the critical coupling strength gc = (U/Jz)c separates a superfluid 
regime Jz ^ U from a Mott insulator regime Jz ^ U (z is the lattice coordination 
number). 

As indicated above, here in contrast we are interested in the competition of Hamil- 
tonian vs. dissipative dynamics. As indicated above, the hopping 7 is a compatible 
energy scale, in the sense that a purely kinetic Hamiltonian has the dissipatively 
targeted |BEC) as an eigenstate. On the other hand, the onsite interaction U coun- 
teracts the off-diagonal order and thus leads to a competition with dissipation of 
strength k. This provides a nonequilibrium analog to the generic purely Hamilto- 
nian equilibrium scenario, in which g = U/kz plays the role of a competition pa- 
rameter - a dominant dissipation g «; 1 supports a condensed steady state, whereas 
dominant interaction g » 1 results in a diagonal density matrix. 

A yet different kind of dynamical phase transitions, which result from the compe- 
tition between different terms of the dissipative Liouvillian, have been anticipated 
in Verstraete et al. (2009), and discussed in more detail in Eisert and Prosen (2010) 
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and Hoening et al. (2010), where in particular the key aspect of criticality in terms 
of diverging length and time scales has been established. Furthermore, our scenario 
is in a sense dual to the dissipative quantum phase transition of a single particle 
on a lattice coupled to a long wavelength heat bath, known to undergo a transition 
from diffusive to localized behavior upon increasing dissipation strength (Schmid, 
1983; Chakravarty et al., 1986; Kampf and Schon, 1987; Chakravarty et al., 1987). 

Theoretical approach - The absence of standard concepts for thermodynamic equi- 
librium, such as the existence of a free energy and associated variational principles, 
makes it necessary to argue directly on the level of the equation of motion (EOM) 
for the density operator, resp. on the associated full set of correlation functions. This 
is in general a formidable task, even numerically intractable in the thermodynamic 
limit in which we are here interested. For this reason, we have developed a gen- 
eralized Gutzwiller mean field approximation scheme, which captures the physics 
in the two well-understood limiting cases g ^ 1,^ » 1, and otherwise provides 
an interpolation scheme. It is implemented by a product ansatz p = ^^pe for 
the full density matrix, such that the reduced local density operators pc = Tr^^p 
are obtained by tracing out all but the fth site. Compared to the standard bosonic 
Gutzwiller procedure for the Bose-Hubbard model at zero temperature, where the 
factorization is implemented for the wave function, it allows for the description of 
mixed state density matrices. It treats the onsite physics exactly, and drops the (con- 
nected) spatial correlations, such that it can be expected to be valid in sufficiently 
high dimensions. The equation of motion for the reduced density operator reads 

drPc = -i[hf,pc\ + £c[pt] , (15) 

where the local mean field Hamiltonian and Liouvillian are given by 



ht = -jYj^{be')b\ + {bl)h)-^ht + "^-Uhtiht - 1), 

(t'\i) 

4 



{t'\l) r,s=\ 

(16) 

hi is in accord with the form of the standard Gutzwiller approach. The addition of 
the chemical potential \x to the Hamiltonian h^ does not change the dynamics, be- 
cause the model conserves the average particle filling n = Yi{{h()IM'^. The freedom 
to fix the chemical potential is necessary to solve the equation dtPss = for the 
steady state of the system (Diehl et al., 2010b; Tomadin et al., 2011). The Liou- 
villian is constructed with the operator valued vector Af = (\,bl,b[,hi), and the 
correlation matrix F^'/ = crVTr^Af ""^''"Af "'"\ with cr = (-1, -1, 1, 1). Note that the 
correlation matrix is p-dependent - this makes the mean field master equation effec- 
tively nonlinear in p. Such a feature is well-known in mean field approximations. 
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e.g. in the Gross-Pitaevski equation, where an A^-body quantum-mechanical linear 
Schrodinger equation is approximated by a non-linear classical field equation. 

The information encoded in Eq. (15) can equivalently be stored in the full set of cor- 
relation functions, resulting in an a priori infinite hierarchy of nonlinear coupled 
equations of motion for the set spanned by the normal ordered expressions {b'^bf) 
for n,m eN and all lattice sites £. This formulation is advantageous in the low den- 
sity limit n «; 1, where we have identified a power counting showing that a closed 
(nonlinear) subset of six correlation functions (i/ff = {be), [b]), {b\b^), c.c), decou- 
ples from the infinite hierarchy. For technical reasons, it is sometimes favorable 
to study the equivalent set of seven connected correlation functions, {i^(, {6b\6bi), 
{5h]), {6b\6b]), c.c), where db^ = be- ifje- This allows to obtain a number of results 
analytically in this limit, such as the condensate fraction as a function of interaction 
strength in the homogeneous limit, and the complete shape of the phase diagram. 

Basic picture for the dynamical quantum phase transition - To better understand 
the phase transition, we consider the limiting cases of vanishing and dominant in- 
teraction. For U = 0, the spontaneous breaking of the U{1) phase symmetry is re- 
flected by an exact steady state solution in terms of a homogeneous coherent state 
pf = |»P)f<^|, with 1^)^ = exp(-n/2)Y,,„[{ne'Y'^)/y[^.]\m)e for any £, together 
with the choice fi = -Jz. The efl'ect of a finite interaction U is best understood us- 
ing a rotating frame transformation on Eq. (15), ViU) = exp[iUhe(he - l)t]. While 
the interaction term is then removed from the Heisenberg commutator, the annihi- 
lation operators become VbeV~^ = 2m Q^piimUt)\m)e{m\be. U therefore rotates the 
phase of each Fock states difl^erently, thus dephasing the coherent state p'fK In con- 
sequence, off-diagonal order will be completely suppressed for sufficiently large U 
and the density matrix takes a diagonal form. Under the assumption of diagonality, 
the master equation reduces to a rate equation 



d,p, = K[(n + \){2bePeb] - {b]be,pe]) (17) 

+n{lb]pebe-{beb],pe})l 

This is the equation for bosons coupled to a thermal reservoir with thermal occu- 
pation n, with thermal state solution p^^'?^ ^ = n"^ /{n + 1)"''^'(5,„,^, where m,k are the 
Fock space indices of the ^th site. At this point two comments are in order. First, 
in contrast to the standard case of an external heat reservoir, the terms n, n -I- 1 are 
intrinsic quantities, meaning that the strongly interacting system provides its own 
effective heat bath. Second, from the solution we note the absence of any distinct 
commensurability effects for integer particle number densities, contrasting the Mott 
scenario at zero temperature. This can be traced back to the fact that in the latter 
case, the suppression of off-diagonal order is additionally constrained by the purity 
of the state, such that (at least on the mean field level) the diagonal pure Mott state 
is the only possible choice. The driven-dissipative system has no such constraint on 
the purity of the state. 
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Fig. 17. Dynamical phase transition: Relaxation dynamics of the condensate fraction to 
the stationary state from an initial fully condensed state as a function of the interaction 
strength U, ioT J = 1.5 K,n = I. Each line corresponds to a stroboscopic snapshot. Inset: 
Near critical evolution reflected by the time-evolution of the logarithmic derivative of the 
order parameter i//{t), for J = 0, n = I, and U < Uc- The early exponential decay (tilted 
crosses) of the initial fully-condensed state is followed by a scaling regime (empty circles) 
with exponent a - 0.5. The final exponential runaway (vertical crosses) indicates a small 
deviation from the critical point. Figure reprinted with permission from Diehl et al. (2010b). 
Copyright 2010 by MacMillan. 

3.2.2 Critical Behavior in Time 

Fig. 17 shows stroboscopically the approach to the steady state in the homogeneous 
limit as a function of interaction strength. In particular, we note the expression of a 
non-analyticity as r ^ oo, characteristic of a second order phase transition. In the 
low density limit, the steady state condensate fraction can be obtained analytically 
and reads 



l-AoP 



1- 



2u'[\+{j + uf) 



l+u^ + j(Su + 6j(l + 2m2) + 24 fu + Sf) 



(18) 



with dimensionless variables u = U/(4kz), j = JI{4k). The boundary between the 
thermal and the condensed phase with varying J,n is shown in Fig. 18 with solid 
lines. 

On general grounds, one expects a critical slowing down at the phase transition 
point when approaching it in time at the critical interaction strength. More pre- 
cisely, the order parameter evolution of the generic form \if/\ ~ exp(-m^t)/f should 
have a vanishing mass or gap term m^ (real part of the lowest eigenvalue), leading to 
a polynomial evolution. The associated scaling of the order parameter is reflected in 
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the plateau regime in the inset of Fig. 17, which sets in after an initial transient and 
is followed by an exponential runaway for a slight deviation from the exact critical 
point. In the low density limit, it is possible to extract the associated dynamical crit- 
ical exponent: At criticality, the order parameter evolution is seen to be governed 
by a cubic dissipative nonlinearity ~ |(/rp, implying solutions \\p\ ^ 1/(4 V^) with 
exponent a = 1/2. This is a mean field result and not indicative of the precise uni- 
versality class of the system, governed by anomalous critical exponents. This issue 
is currently under investigation in a Keldysh path integral approach. Nevertheless, 
already the above result highlights that in our dynamical system, criticality could 
be monitored directly as a function of time, e.g. by stroboscopically measuring the 
condensate fraction. 



3.2.3 Dynamical Instability and Spontaneous Translation Symmetry Breaking 

An intriguing feature of the non-equilibrium stationary state phase diagram is an 
extended region in parameter space, where both the symmetries of phase rotations 
and translations are broken spontaneously, in this sense defining a supersolid phase. 
This state is characterized by a density modulation which is incommensurate with 
the lattice spacing. As illustrated in the phase diagram Fig. 18, the effect occurs 
universally in all density regimes. The plausibility for such a new qualitative effect 
can be understood from the fact that the (bare) dissipation rate /Cq ~ q^ (see below), 
vanishes in the vicinity of the dark state at q = 0: In consequence, there will al- 
ways be a momentum scale where even an arbitrarily weak interaction energy Un 
becomes comparable. In the low density limit, it is possible to describe the phe- 
nomenon analytically, in this way getting insights into the origin of the additional 
phase with translation symmetry breaking. To this end, we work with the closed 
subset of seven correlation functions defined above, which however are time and 
space dependent. Working in a linear response strategy, we linearize around the ho- 
mogeneous steady state solution to study its stability. Upon Fourier transform, we 
obtain a 7 x 7 matrix evolution equation. 

We linearize in time the EOM of Eq. (15), writing the generic connected corre- 
lation function as {Oi;){t) = (<9^)o + 6{0(){t), where {0()q is evaluated on the 
homogeneous steady state of the system. The EOM for the time and space de- 
pendent fluctuations 6{0(){t) = 5Q>({t) is then Fourier transformed, resulting in a 
7x7 matrix evolution equation dt5Q>(^ = Mq^Oq for the correlation functions Oq = 
({6b\, {6b'^)-q, {6b'^6b\, {6b\, {6b^^)-^, {6b^6b\, {6b'^^6b)^^) (We note that the 
fluctuation 6{5b)n {5{6b^)(^ coincides with the fluctuation of the order parameter 
6^|/(^ ((5(Alq), since the average of JZ?, on the initial state vanishes by construction.): 
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Fig. 18. (Color online) Stationary state phase diagramfordifFerentregimes of density: « = 1 
(black), « = 0.1 (red: analytical low density limit calculation; blue: numerical low density 
calculation). The coincidence of analytical and numerical results is enhanced as « ^ 0. All 
regimes of density exhibit the same qualitative features with the three phases discussed in 
the text. Figure adapted from Diehl et al. (2010b). 

where we have separated a slowly evolving sector describing the single particle 
fluctuations and containing the dark state S^i^q = (6tffq, Si/r* ), and a sector *P2,q = 
({6b'^6b)^, {6b^)q, {6b^^)-^, {6b^5b^)(^, {5b^^6b)-^), whose evolution is seen to be lower 
bounded by the scale kh. This matrix is easily diagonalized numerically, with the re- 
sult for the imaginary part of the different branches, describing the damping, shown 
in Fig. 19a. A separation of scales for the lower branches *Fi and the higher ones 
^2 is clearly visible for low momenta q ^ 0, suggesting to integrate out the fast 
modes by adiabatic elimination dtS^j = 0. This results in a renormalization of the 
single particle complex excitation spectrum via the terms involving fractions, 



dt 
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q/ 



where 6q = 7q^ is the kinetic energy and A-q = 2{2n + 1)ko^ the bare dissipative 
spectrum for low momenta. The low-momentum spectrum of this matrix reads 



rq ^ iciql + /Cq, c = A,^LW^^9LW(2z)], 



(21) 



with c the speed of sound. The quadratic q-dependence present without renormal- 
ization correction is modified by a nonanalytic linear contribution, which dominates 
at small momenta and reproduces the shape of the unstable modes obtained via di- 
agonalization in Fig. 19. For a hopping amplitude smaller than the critical value 
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Fig. 19. (Color online) Dynamical instability: a) Damping spectrum as a function of quasi- 
momentum from linear response around a homogeneous state. There are rapidly damping 
branches evolving at 0(K),0(Kn), as well as two slowly evolving branches associated to 
single particle excitation damping. Around the dark state at ^ = 0, a continuum of unstable 
modes appears, b) Numerical evolution of the nonlinear system of correlation functions 
in the low density limit for 800 lattice sites. The color code represents the density profile, 
demonstrating an incommensurate charge density wave stationary phase with characteristic 
wavelength /IcDW- Figure adapted from Diehl et al. (2010b) and Tomadin et al. (2011). 

7c = 9Un/(2z), the speed of sound becomes imaginary, rendering the system un- 
stable. The linear slope of the stability border for small J and U is clearly visible 
from the numerical results in Fig. 18. 

Beyond the unstable point, the linearization strategy around the homogeneous state 
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fails in describing the true steady state of the system. In order to extract the correct 
stationary state in this regime, we resort to a numerical treatment of a large system 
in the low density limit, where the nonlinearities are fully taken into account. The 
result is displayed in Fig. 19b, revealing that the stationary state exhibits charge 
density wave order with characteristic wavelength /Icdw which is set by the in- 
verse of most unstable momentum mode. Generically, it is incommensurate with 
the lattice spacing. The scale characterizing the instability is thus transmuted into 
a physical length scale. The phenomenon is found universally for different system 
sizes, ruling out the possibility of a mere finite size effect. 

At this point, three comments are in order. First, the subtle renormalization effect 
is not captured by a Gross-Pitaevski type approximation scheme and relies on a 
suitable treatment of the higher order correlation functions. Second, the new phase 
emerges at weak coupling already, and for small enough J the homogeneous dis- 
sipative condensate is unstable towards the pattern formation at arbitrarily weak 
interaction. In this weak coupling regime, our approximation scheme is very well 
controlled. Third, the effect relies on the existence of a continuum of modes, and 
thus has a truly many-body origin. In summary, the phase with simultaneous spon- 
taneous breaking of phase rotation and lattice translation symmetry is understood as 
a fluctuation induced beyond (standard) mean field many -body phenomenon, which 
seems quite unique to the dissipative setting. The full phase diagram discussed here 
is shown in Fig. 18. 



3.3 Dissipative D-Wave Paired States for Fermi-Hubbard Quantum Simulation 



Motivation - One of the big experimental challenges in the field of cold atoms is 
the quantum simulation of the ground state of the Fermi-Hubbard model (FHM) 
describing two-component fermions interacting locally and repulsively on the lat- 
tice, whose filling is controlled by a chemical potential. The particular interest in 
this model roots in the fact that it is believed to be a minimal model for the de- 
scription of cuprate high-temperature superconductors. The model has challenged 
theorists for almost thirty years by now, and has proven to be hard to analyze with 
both advanced analytical approaches and numerical techniques. In particular, from 
the theory point of view, so far the d-wave ordered nature of the ground state away 
from half filling, which is observed experimentally, has only the status of a conjec- 
ture. Together with the uncertainty whether the model actually faithfully captures 
the microscopic physics of the cuprates, this situation calls for a quantum simula- 
tion of the FHM ground state in a cold atom context, taking advantage of precise 
microscopic control in such systems. 

This goal still remains very challenging, due to tough requirements on the temper- 
ature in these systems. In fact, the d-wave gap in the cuprates, setting the tempera- 
ture scale to be reached, is only ~ 0.0 ITp (T^ the Fermi temperature), and therefore 



49 



still more than an order of magnitude away from what can currently be reached in 
the lab. Despite impressive progress in this direction (Hofstetter et al., 2002; Kohl 
et al., 2005; Chin et al., 2006; Jordens et al., 2008; Schneider et al., 2008; Jordens 
et al., 2010; Esslinger, 2010), where quantum degeneracy is reached on the lattice, 
new cooling strategies are needed to achieve this goal. The roadmap using dissipa- 
tion state engineering is the following: (i) We dissipatively produce a low entropy 
state that is "close" (in a sense specified below) to the expected ground state of the 
Fermi-Hubbard model away from half filling, (ii) We then construct a suitable adi- 
abatic passage, that consists in slowly switching off the Liouville dynamics while 
ramping up the Hubbard Hamiltonian. 

Here, we will present a mechanism which allows to engineer fermionic paired states 
of arbitrary symmetry, exemplified here for the case of d-wave symmetry, which is 
based on dissipative dynamics alone and works in the absence of any attractive 
conservative forces. The mechanism is based on an interplay of the above mech- 
anism of quasilocal phase locking, and Pauli blocking, thus crucially relying on 
Fermi statistics. A suitable mean field theory, valid for the long-time evolution, has 
a natural interpretation in terms of damping of fermionic quasiparticles and sim- 
plifies the microscopically quartic (interacting) Liouville operator into a quadratic 
one. We then discuss possible implementations and present numerical results for a 
suitable adiabatic passage. 

The state to be prepared - We target BCS-type states, which represent the con- 
ceptually simplest many-body wave functions describing a condensate of A^ paired 
spin- 1/2 fermionic particles. Working on a bipartite square lattice, and assuming 
singlet pairs with zero center-of-mass momentum, we have 

IBCSw) ~ {d^f'^\vsic), (22) 

q iJ 

where Cq o- {c]^) denotes the creation operator for fermions with quasimomentum q 
(on lattice site z) and spin a =T, i, and ^q {(fij) the momentum (relative position) 
wave function of the pairs. We now specialize to a state close to the conjectured 
FHM ground state, in what concerns (i) the symmetries and (ii) the ground state 
energy. For the above pair creation operator d\ the pair wave function 



^q = COS q^ - COS qy or (pij = | ^ pA(5ij+e, + Sij.^J (23) 

A=x,y 

withp_r = -py = 1 ensures the symmetry properties of pairing in the singlet channel 
and the d-wave transformation law ifq^^q^. = -(f-^^.^q^ = <P-q^-q,. under spatial rota- 
tions. The wave function corresponds to the limit of small pairs (see Fig. 20a), and 
phase coherence is granted by the delocalization of these molecular objects. Pairs 
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with such a short internal coherence length appear in the cuprates in the regime 
where strong correlations set in upon approaching half filling. No quantitative state- 
ment can, of course, be made on the energetic proximity of this wavefunction to the 
true FHM ground state. However, the fact that the pairing occurs off-site avoids ex- 
cessive double occupancy (which is energetically unfavorable for the strong repul- 
sive onsite interactions), and makes this state an interesting candidate for quantum 
simulation. 



3. 3. 1 Dissipative Pairing Mechanism 

We now construct a parent Liouvillian, which has the above d-wave state |d) as a 
dark state. In other words, we will construct a set of (non-hermitian) jump operators 
with the property /^'''IBCSat) = 0, where i = I, ...,M (a = x,y,z) represents a posi- 
tion (spin) index (M is the number of sites in the lattice). Due to the product form 
of the dark state wavefunction, a key sufficient condition to fulfill this task is to find 
a set of normal ordered jump operators Jf, which commute with the generator of 
the dark state. 



[Jf,d'^]=Oyi,a. (24) 

The appearance of both indices reflects the need to fix the properties of the state in 
both position and spin space. From a practical point of view, we require the jump 
operators to be quasilocal, number conserving (i.e. [J", N] = 0) and to act on single 
particles only, restricting their class to quasilocal phase rotation invariant fermion 
bilinears. The above condition is very general and thus applicable to wider classes 
of paired, or even more generally, product states. One example discussed in the next 
section is p-wave paired states for spinless fermions. 

We now turn to the construction of the Lindblad operators for the d-wave BCS 
state. To this end, we follow the physical picture that d-wave superconductivity (or 
superfluidity) on a lattice can be viewed as delocalized antiferromagnetic order, ob- 
tained when moving away from half filling (Anderson, 1987; Zhang et al., 1988; 
Gros, 1988; Altman and Auerbach, 2002; Paramekanti et al., 2004). Therefore, 
we will first construct the parent Liouvillian for a Neel state at half filling, which is 
the conceptually simplest (product) wavefunction representing antiferromagnetism, 
and then generalize to the BCS state. There are two Neel states at half filling, re- 
lated by a global spin flip, \N+) = O/eA cj,^^ .,.cjjvac), |N-) = n,e/icj,^ ^cj^|vac) 
with A a sublattice in a two-dimensional bipartite lattice. For later convenience we 
introduce "Neel unit cell operators" S'^^ = cj^^ cr"c] (a = +,ev = {±6;^, ±6^}, and 
two-component spinor c, = (c;,^,c/,j,)), such that the state can be written in eight 
different forms, |N+) = UieASfJyac) = (-1)'^'^ Y\ieBSJ_ Jya.c). We then see that 
the Lindblad operators must obey [j"y,S^j] = for all i,j located on the same 
sublattice A or B, which holds for the set 
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Fig. 20. (Color online) D-wave state and action of the jump operators, a) Symmetries of 
the state: an offsite fermion singlet pair exhibits a characteristic sign change under spatial 
rotations, and is delocalized over the whole lattice. b,c) The dissipative pairing mechanism 
combines (b) Pauli blocking and (c) delocalization via phase locking, b) The action of 
Lindblad operators using Pauli blocking for a Neel state (see text), c) The d-wave state 
results as a delocalization of these pairs away from half filling (shown is a cut along one 
lattice axis). Figure adapted from Diehl et al. (2010c). 



Jty = cLyci,i^AoTB. 



(25) 



The presence of fermionic statistics is essential for the action of the operators j"^, 
as illustrated in Fig. 20b: they generate spin flipping transport according to e.g. 
j^^ = cj^g -,c,-4, not possible when the Neel order is already present. It is then easy 
to prove the uniqueness of the Neel steady state up to double degeneracy: The 
steady state must fulfill the quasilocal condition that for any site occupied by a 
certain spin, its neighboring sites must be filled by opposite spins. For half filling, 
the only states with this property are |N+). The residual twofold degeneracy can be 
lifted by adding a single operator j. 



cT^g (1 + (T^)Ci on an arbitrary site /. 



The Lindblad operators for the d-wave BCS state can now be constructed along a 
similar strategy. First we rewrite the d-wave generator in terms of antiferromagnetic 
unit cell operators 5", 



d' = h Z^^L. - cLyci = f J] D^ D': = j]pji 



(26) 



where p+x = -p±y = 1, and the quasilocal d-wave pair D^ may be seen as the "d- 
wave unit cell operators". This form makes the picture of d-wave superconductivity 
as delocalized antiferromagnetic order transparent, and we note the freedom a = + 
in writing the state. The condition {J'l , Y^j ^^] = (a = (a, z) or (jc, y, z)) is fulfilled 
by 



= E 



PvJiv: 



i" = c^ cr"c- 



(27) 
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which is our main result. Coherence is created by these operators via phase locking 
between adjacent cloverleaves of sites. 

The uniqueness of this state as a stationary state for the Lindblad operators (27) 
is less obvious then in the antiferromagnetic case and we argue based on symme- 
try arguments. Uniqueness is equivalent to the uniqueness of the ground state of 
the associated hermitian Hamiltonian H = VYji,a=±,z-fT •^'i ^^^ V > Q. The state 
generated by (27) shares the Hamiltonian symmetries of global phase and spin ro- 
tations, and translation invariance. Assuming that no other symmetries exist, we 
then expect the ground state to be unique. The full set [J'^] is necessary for unique- 
ness: Omitting e.g. {J^) generates an additional discrete symmetry in H resulting 
in two-fold ground state degeneracy. We confirmed these results with small scale 
numerical simulations for periodic boundary conditions, cf. Fig. 21. We note that 
the above construction method allows us to find "parent" Lindblad operators for a 
much wider class of BCS-type states, see Yi et al. (201 1). 



3.3.2 Dissipative Gap 

A remarkable feature of the dissipative dynamics defined with the set of operators 
(27) is the emergence of a "dissipative gap" in the long time evolution of the master 
equation. Such a dissipative gap is a minimal damping rate which crucially remains 
finite in the thermodynamic limit. The phenomenon is a dissipative counterpart of 
a coherent gap suppressing single particle fermion excitations in a BCS superfluid, 
where it is a characteristic feature of the low energy effective theory. 

The dissipative gap can be established in a mean field theory which is controlled 
by the proximity to the exactly known stationary dark state. For this purpose it is 
convenient to give up exact particle number conservation and to work with fixed 
phase coherent states IBCSti) = A/'"^^^exp(e'''J"'')|vac) instead of the fixed number 
states IBCSat) (Leggett, 2006), where N = OqCl + V'q) ensures the normalization. 
The equivalence of these approaches in the thermodynamic limit is granted by the 
fact that the relative number fluctuations in BCS coherent states scale ~ 1 / VA^, 
where A'^ is the number of degrees of freedom in the system. The density matrix 
for the coherent states factorizes in momentum space exp(e"'<i'^)|vac) = OqCl + 
^'^V'qfq t'^Li^l^^''^- ^^ ^^^^ times, we can make use of this factorization property 
and expand the state around jBCS^i), implemented with the ansatzp = OqPq^ where 
Pq contains the mode pair +(q, a) necessary to describe pairing. We then find a 
linearized evolution equation for the density operator. 



JJp] = 2 Kq[7q,<Tpyi^ - ^{ri,^rq,^,p}], (28) 

with quasiparticle Lindblad operators and momentum dependent damping rate given 
by 
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Fig. 21. (Color online) Uniqueness of the d-wave dark state for the master equation with 
Lindblad operators from Eq. (27): Fidelity to the d-wave BCS state, (BCSAr|p|BCSiv) for 4 
atoms on a 4x4 grid, showing exponential convergence from a completely mixed state to a 
pure state. Dashed lines denote sampling errors. (Inset): Entropy evolution for four atoms 
on a 4x1 lattice. Figure adapted from Diehl et al. (2010c). 

7q,o- = (1 + fl)'^'^ ic~q,,T + Sa-ifqcl^-a-), 
Kq = Kh(l + V?q) > KH, 



(29) 



with Sf = -l,si = 1, the wavefunction specified in Eq. (23), and the value h = 
2 J T^ iri^ ~ ^'^^ dictated by the presence of nonzero mean fields resulting 
from a coupling to other momentum modes, and the proximity to the final state. 

The linearized Lindblad operators have analogous properties to quasiparticle op- 
erators familiar from interaction pairing problems: (i) They annihilate the (unique) 
steady state yq^o-IBCSg) = 0; (ii) they obey the Dirac algebra {yq,^, y^,^,} = 6q,q>Sa;o-' 
and zero otherwise; and (iii) in consequence are related to the original fermions via 
a canonical transformation. 

Physically, the dissipative gap kh implies an exponential approach to the steady 
d-wave BCS state for long times. This is easily seen in a quantum trajectory rep- 
resentation of the master equation, where the time evolution of the system is de- 
scribed by a stochastic system wavefunction \tffit)) undergoing a time evolution 
with non-hermitian "effective" Hamiltonian \il/(t)) = e~'"^<^^\il/(0)) / ||. . .|| (//eff = 
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H-'iK Yiia Jf^Jf here) punctuated with rate k Il7VI(A(0)lP by quantum jumps |(A(0) -^ 
jV|i/r(0)/ll- • -11 such that p(t) = (|iA(0)(iA(OI)stoch (see, e.g., Gardiner and Zoller 
(1999)). We thus see that (i) the BCS state is a dark state of the dissipative dy- 
namics in that jV|BCSa?) = implies that no quantum jump will ever occur, i.e. the 
state remains in IBCSat), and (ii) states near |BCS/v) decay exponentially with rate 
lower-bounded by the dissipative gap. 

This dissipatively gapped behavior strongly contrasts the bosonic case, where the 
dissipation is gapless as we have seen above, in the sense that /Cq ~ q^ for q ^ 0. 
One crucial difference between the bosonic and fermionic evolutions is then the 
fact that many-body observables involving a continuum of modes behave polyno- 
mially in the boson case, due to the slow decay in the vicinity of the dark state. 
For fermions instead, the dark state property is not encoded in a zero of the decay 
rate, but rather in the annihilation property of the linearized Lindblad operators on 
a nontrivial BCS vacuum. In this case, even many-body observables will relax ex- 
ponentially. More generally, the generation of a finite gap scale at long times makes 
the fermionic dissipatively induced phases potentially more stable than the bosonic 
ones, as one may compare competing energy or rate scales to that finite scale. 

This convergence to a unique pure state is illustrated in Fig. 21 using numerical 
quantum trajectory simulations for small systems. We show the fidelity of the BCS 
state for a small 2D grid as a function of time, computed for the full density matrix 
via the quantum trajectories method. The inset shows the entropy evolution for a 
small ID system (where one direction of the d-wave cloverleaf is simply omitted). 



3.3.3 State Preparation 

Implementation with alkaline earth-like atoms - The conceptually simple quasilo- 
cal and number-conserving form of J" raises the possibility to realize dissipation 
induced pairing via reservoir engineering with cold atoms. We illustrate this in ID, 
taking the example of J^ = (cj^j -. -I- c[ j^)cj^. Implementation requires (i) a spin 
flip, (ii) a spatial redistribution of the atom onto sites neighbouring the central one, 
and (iii) a dissipative process which preserves the coherence over several lattice 
sites. These ingredients can be met using alkaline earth-like atoms (Ye et al., 2008; 
Reichenbach and Deutsch, 2007; Daley et al., 2008; Gorshkov et al., 2010) with 
nuclear spin (e.g., 7=1/2 for "'Yb), and a long-lived metastable ^Pq manifold 
as a physical basis, see Daley (2011) for a recent review. In this setting, one can 
construct a stroboscopic implementation, where the action of each J" is realized 
successively. The level scheme and the spin flip process are described in Fig. 22a. 
There we concentrate on the spatial redistribution of the atoms using the fact that 
the ^Po states can be trapped independently of the ground ^So manifold. The ^Pq 
state is trapped in a lattice of three times the period as that for the 'So state, defin- 
ing blocks of three sites in the original lattice. Using this, any J, atom in 'So on the 
central site is excited to the T state of the ^Pq manifold. By adding an additional 
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Fig. 22. (Color online) a) Level scheme for physical implementation. The spin flip opera- 
tion is implemented via ofF-resonant coherent coupling to the ^Pq manifold with circularly 
polarized light (red arrows). The long lived ^Pq states are coupled to the ^Pi level in a 
two-photon process, from which spontaneous emission into a cavity is induced, leading 
back to the ^^o manifold encoding the physical fermionic states, b) Adiabatic passage 
connecting the mean field d-wave state with the ground state of the FHM on a 2 x 6 lad- 
der with 4 atoms with parent Hamiltonian //p (see text). Evolution of fidelity of the in- 
stantaneous system state with respect to the final ground state of the FHM is calculated, 
(inset): Parameters hopping J, onsite interaction IJ of the Fermi-Hubbard Hamiltonian 
H = -J Z<(,;>,o- (^ia-^jcr + U J^i ^|t^'T^;i ^'i' ^^^ '^^c parent Hamiltonian strength V, as a func- 
tion of time in units of the maximal final interaction strength [/„,. Figure adapted from Yi 
etal. (2011). 

potential, the traps for ''Pq are coherently divided so that atoms confined in them 
overlap the right and left sites of the original block. Decay is induced by coupling 
atoms in the ^Po state ofF-resonantly to the ^Pi state, as depicted in Fig. 22a, with 
coupling strength Q, and detuning A. By coupling the ^Sq-^Pi transition to a cav- 
ity mode with linewidth Y and vacuum Rabi frequency g, the decay is coherent 
over the triple of sites. In the limit A » 11 and F » "£' ^^ effective decay rate 

Teff = -^ ~ 9kHz results for typical parameters. Fermi statistics will be respected 
in this process, as long as the atoms remain in the lowest band. This operation can 
be performed in parallel for different triples, and needs to be repeated with the su- 
perlattice shifted for other central sites. Similar operations combined with rotations 
of the nuclear spin before and after these operations allows implementation of J7 
and J^. In 2D 3x3 plaquettes are defined by the appropriate superlattice potential 
for the ''Po level, and the adiabatic manipulation of the potential has to be adjusted 
to ensure the correct relative phases for atoms transported in orthogonal directions. 
Such a digital or stroboscopic scheme is rather demanding in the context of cold 
atoms, and most of the complication comes from the need to fix the spin quantum 
number. Below, we discuss spinless fermions and see that there, an "analog" im- 
plementation along the lines of Sect. 3.1 with continuous driving and dissipation is 
possible. 
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Adiabatic Passage - To reach the ground state of the FHM in small scale numerical 
simulations, we found it efficient to introduce in addition to the parent Liouvillian 
a parent Hamiltonian Hp = V 2,> J"^Jf, which has the above d-wave state as the 
exact unique (fixed number) ground state for V > 0, and which could be obtained by 
replacing the decay step into the cavity by induced interactions between atoms. The 
result of the numerical calculation is reported in Fig. 22b, where convergence to the 
FHM ground state is clearly seen. In a large system, one should additionally be able 
to take advantage of the fact that (i) in the initial stages the system is protected by a 
gap ~ 0.12V, and (ii) the d-wave state has identical symmetry and similar energy to 
the conjectured Fermi-Hubbard ground state away from half filling. Thus, a d-wave 
superfluid gap protection is present through the whole passage path, since no phase 
transition is crossed. 



3.4 Dissipative Topological States of Fermions 



Motivation - Topological phases of matter exhibit ordering phenomena beyond the 
Landau paradigm, where order is described by local order parameters. Instead, 
these phases are characterized by nonlocal order parameters, the topological in- 
variants (Hasan and Kane, 2010; Qi and Zhang, 201 1). Observable physical mani- 
festations of topological order emerge when these systems are subject to boundary 
conditions in space, such as the appearance of Majorana modes localized to suitably 
designed edges in certain one- or two-dimensional superfluids (Kitaev, 2001; Read 
and Green, 2000). These modes are robust against large classes of environmental 
perturbations and imperfections. This gives them a potentially high practical rele- 
vance, and they are discussed as candidates for providing the building blocks for 
topologically protected quantum memories and computations (Nayak et al., 2008). 

So far, the concept of topological order and its physical consequences have been 
discussed mainly in a Hamiltonian ground state context. Motivated by the prospects 
of combining topological protection with a targeted dissipative engineering of the 
corresponding states, in Diehl et al. (2011); Bardyn et al. (2012) we have shown 
how such concepts and phenomena manifest themselves in systems governed by 
driven-dissipative Lindblad dynamics. Here we will give a brief review of these re- 
sults, focusing on the simplest paradigmatic model discussed in Diehl et al. (201 1), 
a dissipative quantum wire of spinless atomic fermions. This model is the counter- 
part of Kitaev's quantum wire, which provides a minimal one-dimensional model 
for topological order, and hosts Majorana edge modes in a finite wire geometry. In 
particular, we establish dissipative Majorana modes, and discuss their interpreta- 
tion in terms of a nonlocal decoherence free subspace. We give an argument for the 
nonabelian exchange statistics, and sketch the construction of a topological invari- 
ant for density matrices corresponding to mixed states pinpointing the topological 
origin of the edge modes. We also highlight a phase transition induced by "loss of 
topology" which has no Hamiltonian counterpart. Beyond these theoretical find- 
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ings, we argue that due to the spinless nature of the atomic constituents, an im- 
plementation along the lines of Sect. 3.1 is possible. Remarkably, all that needs to 
be done is to replace the bosonic operators in Eq. (11) by spinless fermionic ones, 
and to put proper boundaries using the new experimental tools offered by singe-site 
addressability (Bakr et al., 2010a,b; Sherson et al., 2010; Weitenberg et al., 2011). 
Together with practical preparation protocols and detection schemes (Kraus et al., 
2012), this makes dissipative state engineering an attractive route for realizing Ma- 
jorana physics in the lab. 

Topological quantum wire in Hamiltonian setting - Before embarking the construc- 
tion of a dissipative quantum wire, we first recapitulate briefly Kitaev's Hamiltonian 
scenario. We discuss spinless fermions ai,a] on a finite chain of A^ sites / described 
by a quadratic Hamiltonian // = !], \-Ja\ai^i - Aa,a,+i -I- h.c. j - //a- a,- , with hop- 
ping amplitude 7, a pairing term with order parameter A, and a chemical potential 
jj.. The topologically non-trivial phase of the model is best illustrated for parameters 
7 = |A| and ju = 0, where the Hamiltonian simplifies to 



N-l I , 1 N-\ 

// = 27 ^ \a]ai - 2 = ^-^ S "'^^ ^2'+i • (30) 



Here we write the Hamiltonian in a complex Bogoliubov basis defined with quasilo- 
cal fermionic quasiparticle operators a,, and in terms of Majorana operators c,, 
which are given by the quadrature components of the original complex fermion op- 
erators a,- = ^ {c2i - ic2,-i), respectively. With these preparations, we collect some 
key properties of this model: The bulk properties are most clearly revealed in the 
complex Bogoliubov basis, where the Hamiltonian is diagonal: The ground state is 
determined by the condition a,|G) = for all /, and the bulk describes a fermionic 
BCS-type p-wave superfluid with a bulk spectral gap, which for the above pa- 
rameter choice equals the constant dispersion 6/t = 27. The Majorana represen- 
tation instead gives rise to a picture of the bulk in terms of pairing of Majoranas 
from different physical sites. In view of the edge physics, the absence of the term 
2{a\aN-ll2) = ic2/vCi for a finite wire indicates the existence of a two-dimensional 
zero energy fermionic subspace spanned by |a) e {|0), |1) = aj^|0)}, which is highly 
non-local in terms of the complex fermions. In contrast, in the Majorana basis the 
situation is described in terms of two Majorana edge modes y^ = ci {jr = C2n), 
which are completely localized on the leftmost (rightmost) Majorana site 1 {2N), 
describing "half" a fermion each. These edge modes remain exponentially local- 
ized in the whole parameter regime -27 < p < 27, however leaking more and 
more strongly into the wire when approaching the critical values. Their existence is 
robust against perturbations such as disorder, which can be traced back to the bulk 
gap in connection with their topological origin (Kitaev, 2001). 
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3.4.1 Dissipative Topological Quantum Wire 

i) Bulk properties - In view of constructing an open system analog of the above 
scenario, we consider a purely dissipative {H = 0) Lindblad master equation of the 
form of Eq. (1) for spinless fermions in a chain with A'^ sites and rate k. We choose 
the Lindblad jump operators j, as the above Bogoliubov quasiparticle operators, 
with the explicit form 

ji = ai = -{ai + a]-aM+a]^^), {i = \,. . . ,N - \). (31) 

These Lindblad operators are quasi-local superpositions of annihilation and cre- 
ation operators, leading to a Liouville operator which is quadratic in the fermions, 
and act on the links of each pair of lattice sites (see Fig. 23a). We indicate below 
how such a setting emerges naturally in the long-time evolution of a microscopi- 
cally number conserving (quartic) Liouville dynamics, relying on a mean-field the- 
ory as discussed in Sect. 3.3.2, and taking advantage of the quasilocal nature of the 
target Lindblad operators. Crucially, the ground state condition a/|G) = now plays 
the role of a dark state condition. Since the operators j, obey the Dirac algebra, in 
a translation invariant setting this dark state is unique and pure. In particular, the 
bulk of the system cools to the p-wave superfluid ground state of the Hamiltonian 
(30). The approach to this steady state is governed by the damping spectrum of the 
Liouvillian X. In analogy to the Hamiltonian gap in Kitaev's model, diagonality 
of X in the a, now implies a flat damping spectrum Kk = k, and in particular the 
existence of a dissipative gap. 

ii) Edge modes as nonlocal decoherence free subspace - For a finite wire we 
find dissipative zero modes related to the absence of the Lindblad operator aN. 
More precisely, there exists a subspace spanned by the edge-localized Majorana 
modes a^ = ^ (Itl + Tr), with the above Fock basis lor) e {|0),|1)}, which is 
decoupled from dissipation, i.e. d,Pap{t) = with p^p = {a\p\P). These dissipa- 
tive edge modes are readily revealed in solutions of the master equation defined 
with jump operators Eq. (31). The fact that the master equation is quadratic in 
the fermion operators implies solutions in terms of Gaussian density operators 
p{t) ~ exp -^c^G(f)c . Here we have defined a column vector c of the 2N Ma- 
jorana operators, and G is a real antisymmetric matrix related to the correlation 
matrix Yab{t) = ^{[Ca,Ch]) = i[tanh(iG/2)]afo, which equally is real and antisym- 
metric. Writing the Lindblad operators in the Majorana basis, j, = l]c, j] = c^l*, 
such that the Liouvillian parameters are encoded in a hermitian 2N x 2N matrix 
M = YjiU ® l\, this covariance matrix obeys the dissipation-fluctuation equation 
(Prosen, 2008, 2010; Eisert and Prosen, 2010), d,T = -{X,T} + Y, with real matri- 
ces X = 2ReM = X'^ and Y = 4ImM = -Y^ . Physically, the matrix X describes 
damping, while the matrix Y is related to fluctuations in the stationary state, deter- 
mined by {X, T} = Y. Note that Y corresponds to the first quantized description of 
the efi'ective Hamiltonian associated to the master equation: Due to Fermi statis- 
tics, only the antisymmetric part of M contains nontrivial information, and thus 
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Fig. 23. (Color online) a) The Lindblad operators act on each link of the finite wire, in 
this way isolating the edge mode subspace described by yi, Jr, which together define the 
Hilbert space of one complex fermion (see text). The bulk (blue shaded) is cooled to a 
p-wave superfluid, with pairing links between different physical sites established dissipa- 
tively. b) Visualization of the winding number v for chirally symmetric mixed states, char- 
acterizing the mapping from the Brillouin zone - S ^ to the vector m^, which due to chiral 
symmetry is constrained to a great circle - S^. For pure states, it is furthermore pinned to 
unit length (large circle). Tuning the Liouville parameters can destroy the purity and de- 
forms the circle to an ellipse (blue), while the topological invariant remains well defined. 
A phase transition occurs when the ellipse shrinks to a line (dark line). The values of 6 are: 
7r/4 (large circle); 1.97r/4 (ellipse); njl (line). Figure adapted from Diehl et al. (2011). 

^eff = Yii i]ii = ^c^Yc. Writing F = F -I- ST, the approach to steady state is gov- 
erned by dt6Y = -{X,6Y}, i.e., the eigenvalues of the positive semi-definite matrix 
X give the damping spectrum. The "dark" nonlocal subspace of edge modes, decou- 
pled from dissipation, is thus associated with the subspace of zero eigenvalues of 
the damping matrix X. We refer to Bardyn et al. (2012) for a more comprehensive 
discussion of the roles of X and Y. 

in) Bulk-edge dynamics and dissipative isolation - In a spectral decomposition 
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^ = Hr ^r\r){r\, and identifying by greek subscripts the zero eigenvalues subspace, 
we can write 



dt 



1 ff/? 1 as 



f 



-(TAXs 



(32) 



While the bulk (rs sector) damps out to the steady state by dissipative evolution, the 
density matrix in the edge mode subspace (a/S sector) does not evolve and therefore 
preserves its initial correlations. The coupling density matrix elements (mixed sec- 
tors) damp out according to Trfj(t) = e~'*'Tr^(0). In the presence of a dissipative gap 
as in the example above, this fadeout of correlations is exponentially fast, leading 
to a dynamical decoupling of the edge subspace and the bulk. 

In summary, we arrive at the physical picture that dissipative evolution cools the 
bulk into a p-wave superfluid, and thereby isolates the edge mode subspace, p(t -^ 
oo) -^ Pedge ® Pbuik, providing a highly nonlocal decoherence free subspace (Lidar 
etal., 1998). 

So far, we did not yet address the preparation of the edge mode subspace. Generi- 
cally, when starting from a wire geometry, the initial edge mode subspace is strongly 
mixed. Since its correlations are preserved during dissipative evolution, it thus will 
be useless e.g. as a building block for a qubit (Note that this property is also shared 
with a Hamiltonian setting, where the equilibrium density matrix peq ~ Q-^Kk-aT)^ ^^ 
the Boltzmann constant, is purified by lowering the temperature. The subspace of 
this density matrix associated to the zero modes of H is not purified by such cool- 
ing.). Therefore, in Kraus et al. (2012) we discuss a scheme where the starting point 
is a ring geometry, where the stationary state is unique and has even parity, since 
it corresponds to a paired state of fermions. The ring is then adiabatically "cut" by 
removing dissipative links quasi-locally. In this way, it is possible to obtain a pure 
Majorana subspace with non-local edge-edge correlations. 



3.4.2 Nonabelian Character of Dissipative Majorana Modes 

There is a simple and general argument for the nonabelian exchange statistics of 
dissipative Majorana modes, highlighting the universality of this property that holds 
beyond the Hamiltonian setting. Consider the time evolution of the density ma- 
trix in a co-moving basis \a{t)) = U{t)\a{0)) which follows the decoherence free 
subspace of edge modes, i.e. preserves pap = 0. Demanding normalization of the 
instantaneous basis for all times, {b{t)\a{t)) = 6ah, this yields 

— p = -i[A,p] + ^ \a)pab{b\, (33) 

a,b 
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with the hermitean connection operator A = itj'^U and Pab = {o-{t)\dtp\b{t)) the 
time evolution in the instantaneous basis. The Heisenberg commutator reflects the 
emergence of a gauge structure (Berry, 1984; Simon, 1983; Wilczek and Zee, 1984; 
Pachos et al., 1999; CaroUo et al., 2003) in the density matrix formalism, which 
appears independently of what kind of dynamics - unitary or dissipative - generates 
the physical time evolution, represented by the second contribution to the above 
equation. We note that an adiabaticity condition Qjkq «c 1 on the rate of parameter 
changes vs. the bulk dissipative gap has to be accommodated in order to keep the 
protected subspace. Since the subspace has no intrinsic evolution, this provides 
a natural separation of time scales which prevents the decoherence-free subspace 
from being left, a phenomenon sometimes referred to as the Quantum Zeno effect 
(Beige et al., 2000). 

Starting from this understanding, one can now construct adiabatic local parameter 
changes in the Liouvillian at the edges of a chain to perform elementary dissipa- 
tive Majorana moves. Applying such procedure sequentially, and operating on a 
T-junction in full analogy to the proposal by Alicea et al. (2011) for Hamiltonian 
ground states in order to exchange the two modes while permanently keeping them 
sufliciently far apart from each other, the unitary braiding matrix describing the 
process is 5,^ = exp(|7;7yj for two Majorana modes i, j. This demonstrates non- 
abelian statistics since [5,^, Bj^] + for i + j. Here we use that the above general 
considerations are not restricted to a single quantum wire but apply to more general 
quantum wire networks. 



3.4.3 Topological Order in Density Matrices 

Density matrix topological invariant - In numerical calculations we have verified 
that the Majorana modes are robust under wide classes of translation-invariance 
breaking perturbations such as random local variations of the Lindblad operators 
of Eq. (31), suggesting a topological origin. Indeed, we can connect the existence 
of the edge modes to topological order in the bulk of the stationary state. This is 
achieved by constructing a topological invariant for the distinction of topologically 
inequivalent states. This classification is formulated in terms of the density matrix 
alone and does not rely on the existence of a Hamiltonian or on the purity of the 
state, in contrast to existing constructions. 

As shown in Diehl et al. (201 1), the topological information of the stationary state 
of a Gaussian translationally invariant Liouville evolution is encoded in the even 
occupation subspace of each momentum mode pair +k, p2k °^ ^(1 + n^^cr), where a 
is the vector of Pauli matrices and fik is a real three-component vector < In^l < 1. 
The special case of pure states corresponds to pi = pk, i.e. |n^| = 1 for all k. In 
the more general case, once the vector n^t is nonzero for all k, a normalized vector 
^k = \nkV^nk can be introduced. This then defines a mapping S^ -^ S^, where S^ 
is the circle defined with the Brillouin zone -n < k < n with identified end points 
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k = +71 as usual, and the unit sphere 5^ is given by the end points of n^, as illus- 
trated in Fig. 23. This mapping, however, is generically topologically trivial, with 
corresponding homotopy group ni(S^) = 0, since a circle can always be continu- 
ously shrunk into a point on the sphere. In order to introduce a nontrivial topology, 
we therefore need an additional constraint on n^ . In our setting, motivated by Ki- 
taev's model Hamiltonian (Kitaev, 2001), this is provided by the chiral symmetry 
(Altland and Zimbauer, 1997; Ryu et al., 2010). In terms of the density matrix, 
the latter is equivalent to the existence of a ^-independent unitary matrix S with 
S^ = 1, which anticommutes with the traceless part of the density matrix {fikcr in 
our case): l^ntcrll = -Hkcr. This condition can be turned into a geometric one, by 
representing the matrix Z in terms of a constant unit vector a, Z = aa. The chiral 
symmetry condition then translates into an orthogonality condition /?/.« = for all 
k. The end point of n^ is now pinned to a great circle S ' on the sphere such that the 
vector n^ defines a mapping S^ ^ S^ from the Brillouin zone into a circle, see Fig. 
23b. The corresponding homotopy group is now nontrivial, ni{S^) = Z, and such 
mappings are divided into different topological classes distinguished by an integer 
topological invariant, the winding number, with the explicit form 

V = — I dkd ■ (uj^x dA) e Z. (34) 

Geometrically, v counts the number of times the unit vector % winds around the 
origin when k goes across the Brillouin zone. Crucially, the resulting topological 
distinction of diff"erent density matrices for translationally invariant, chirally sym- 
metric Gaussian systems works without restriction on the purity of the state. Using 
bulk-edge correspondence established for Hamiltonian settings (Hatsugai, 1993; 
Kitaev, 2006), a nonzero value of the invariant would imply the existence of edge 
modes as found above. However, in a general dissipative setting it is possible to 
break this bulk-edge correspondence. For a discussion of this subject, and interest- 
ing consequences of it, we refer to Bardyn et al. (2012). 

Phase transition by "loss of topology" - In Fig. 23b we illustrate a situation de- 
scribed by a one-parameter deformation of the vector n^{6). This is induced by a 
corresponding deformation on the Lindblad operators according to ji{9) = 4= (sin 6 (a.- 

aj+i) + cos 6 (ai + a]_^^)), where Eq. (31) is reproduced for 9 = n/4. For this deforma- 
tion, the purity is not conserved while preserving the chiral symmetry, reflected in 
the fact that the vector in general lies on an ellipsis < In^:] < 1 . Topological order 
is meaningfully defined as long as the first inequality is strict as discussed above, 
i.e. as long as there is a "purity gap". However, at the points 6 = 9^ = ns/2 {s 
integer), not only the direction of a but also the topological invariant is not defined, 
since fik, aligned in the y-direction for all k, has zeroes and the purity gap closes: 
nj.^o,;r = 0, meaning physically that these modes are in a completely mixed state. 
The "loss" of topology at 9 = 9^ can be viewed as a non-equilibrium topological 
phase transition(Rudner and Levitov, 2009; Lindner et al., 2011; Kitagawa et al., 
2010) as a result of changing the Liouville parameters: 9 = 9^ also implies a clos- 
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ing of the dissipative gap in the damping spectrum, which leads to critical behavior 
manifesting itself via diverging time scales, resulting e.g in an algebraic approach 
to steady state (as opposed to exponential behavior away from criticality) (Diehl 
et al., 2008; Kraus et al., 2008; Verstraete et al., 2009; Diehl et al., 2010b; Eisert 
and Prosen, 2010). We emphasize that the symmetry pattern of the steady state is 
identical on both sides of the transition, ruling out a conventional Landau-Ginzburg 
type transition and underpinning the topological nature of the transition. 



3.4.4 Physical Implementation 

As mentioned above, a physical implementation of this scenario is provided by a 
microscopically number conserving Liouville dynamics as discussed in Sect. 3.1.1, 
with jump operators of the form Eq. (11), where boson operators are replaced by 
spinless fermionic ones. We note that in this implementation setting, the role of 
the bath is played by the bosonic atoms from a surrounding BEC, and originates 
microscopically from standard contact density-density interactions, thus imposing 
a natural parity conservation for the fermionic system constituents due to fermionic 
superselection rules. This contrasts potential solid state realizations, where the en- 
vironmental degrees of freedom are fermionic as well. Explicitly, we choose 

Ji = -(a] + aj^j)(a,- - a,+i) = CU,-. (35) 

From a formal point of view, the sequence of annihilation (A, = |(a, - a;+i)) and 
creation (Cj = ^(aj+a]_^i)) part, gives rise to dissipative pairing of spinless fermions 
in the absence of any conservative forces, in complete analogy to the discussion 
for the spinful case in Sect. 3.3.1. In the present case, the mean field construction 
outlined above can be simplified. It can be shown (Diehl et al., 2011) that in the 
long-time and thermodynamic limit, the following general relation between fixed 
number (7,) and fixed phase (j,) Lindblad operators holds, 

Ji = CjAi^ji = Cj+U (36) 

The relation to the Majorana operators is now apparent: It is precisely Kitaev's 
quasiparticle operators which are obtained as effective Lindblad operators in the 
late time evolution, j, = a,. The role of phase fluctuations remains to be investi- 
gated. The explicit mean field calculation shows that a master equation with jump 
operators Eq. (31) is produced, with effective dissipative rate k = k/S. An analy- 
sis of the leading imperfections shows that they preserve the chiral symmetry, and 
so keep the system in the above described topological class. We furthermore em- 
phasize that the practical simplifications in view of engineering such dissipative 
dynamics in the lab compared to the stroboscopic implementation of Sect. 3.3.3 is 
mainly due to the fact that the spin quantum numbers do not have to be fixed in the 
present case. 
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Recently, we have also investigated two-dimensional dissipatively induced fermionic 
paired states with p^ + ipy order parameter (Bardyn et al., 2012). Intriguingly, in 
such systems we established a mechanism that guarantees the existence of a single 
localized Majorana mode at the core of a dissipative vortex in a phase with vanish- 
ing bulk topological invariant. This phenomenon ultimately relies on a violation of 
the bulk-edge correspondence which is unique to the dissipative dynamics and has 
no Hamiltonian counterpart. The Majorana modes could be generated dynamically 
with the implementation strategy outlined here by additionally imprinting optical 
angular momentum onto the matter system (Brachmann et al., 2011), potentially 
circumventing the need of single-site addressability. 



4 Outlook 



In the present work we have summarized recent advances in digital quantum simu- 
lation and engineering of open many-body systems with atoms and ions, where our 
main emphasis has been on presenting new concepts and tools. We conclude our 
discussion with few remarks on open theoretical and experimental problems and 
challenges. 

With regard to the digital quantum simulation approach discussed in Section 2, 
the described experiments realized with trapped ions (Lanyon et al., 2011; Bar- 
reiro et al., 2011) demonstrate in principle the feasibility of the digital simulation 
approach for the study of open many -particle quantum systems. They have been 
carried out in setups of linear ion chains and are, in their present form, not im- 
mediately scalable to large systems. However, similar protocols can be realized 
in scalable and two-dimensional ion-trap architectures, whose development is cur- 
rently at the center of an intense research effort (Blakestad et al., 2009; Home et al., 
2009; Hensinger et al., 2006; Schmied et al., 2009; Clark et al., 2009). 

In view of the big challenge of scaling up the simulations to larger systems, the 
Rydberg-based simulator architecture with cold atoms in optical lattices (Weimer 
et al., 2010) provides an a priori scalable simulation platform. Especially in view 
of the recent experimental achievement of the first entangling Rydberg gates (Isen- 
hower et al., 2010; WiUc et al., 2010) and single-site addressability (Bakr et al., 
2010a; Sherson et al., 2010; Bakr et al., 2010b; Weitenberg et al., 201 1), it seems to 
be a promising route towards large-system digital quantum simulators with control 
over some tens to hundred qubits (spins). This would outperform state-of-the-art 
classical numerical simulation techniques. However, it remains to be seen if neutral 
atoms or other competing simulation platforms will be able to achieve the remark- 
able fidelities of quantum gate operations demonstrated with ions (Lanyon et al., 
2011; Barreiro et al., 2011). In fact, the concepts discussed here for trapped ions 
and Rydberg atoms can be readily adapted to other physical simulation platforms 
ranging from optical, atomic and molecular systems to solid-state devices (Ladd 
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et al., 2010; O'Brien, 2007; Clarke and Wilhelm, 2008; Wrachtrup and Jelezko, 
2006; Hanson et al., 2007; Vandersypen and Chuang, 2004). 

From a fundamental point of view, it will be most interesting to connect the driven- 
dissipative ensembles discussed in Section 3 to other fields, such as non-equilibrium 
statistical mechanics. For example, as known from classical problems in this con- 
text (see, e.g., Altland and Simons (2010)), strong non-equilibrium drive can give 
rise to new universality classes beyond those known in thermodynamic equilib- 
rium. It seems plausible that similar phenomena could be present in our systems 
as well, possibly enriched by their intrinsic quantum mechanical character. More 
broadly speaking, the goal is the identification of universal hallmark signatures for 
the intrinsic non-equilibrium nature of these systems. 

From a practical perspective, the recent experiments with atomic ensembles (Krauter 
et al., 2011) and trapped ions (Barreiro et al., 2011; Lanyon et al., 2011) suggest 
that a strong dissipative drive can protect against ubiquitous unwanted decoher- 
ence mechanisms - while a system with dominant unitary dynamics alone is sen- 
sitive to decoherence. This sparks the more general question if systematic criteria 
for the stability of many-body states under competing unitary and dissipative dy- 
namics can be established, starting from the promising results on the existence 
of dissipative gaps described above. A general scenario of dissipative protection 
clearly would give a high practical relevance to dissipative quantum state engineer- 
ing. Ultimately, if these questions can be answered positively, it will be intriguing 
to investigate whether the robustness benefits of dissipative quantum computation 
(Verstraete et al., 2009) and memories (Pastawski et al., 201 1), as well as topolog- 
ical quantum computation (Nayak et al., 2008), can be sensibly combined in one 
unified setting. Clearly, answering such questions also requires the development of 
new theoretical tools. A promising avenue is provided by a Keldysh functional inte- 
gral approach (Kamenev and Levchenko, 2009), within which the powerful toolbox 
of advanced field theoretical methods could be leveraged over to driven-dissipative 
many-body systems. 

Furthermore, in view of quantum engineering, it is an important goal to extend 
the scope of many-body physics with driven-dissipative ensembles to new phys- 
ical platforms. This concerns not only trapped ion systems, but also microcavity 
arrays, which have a strong potential of being developed into true many-body sce- 
narios in the future (Mariantoni et al., 2011). Each of these systems will also add 
new theoretical challenges, such as the intrinsic non-number conserving nature of 
systems whose basic constituents are photons. In addition, it will be intriguing to 
explore the theoretical crosshnks between analog and digital quantum simulations 
in a many-body context. 

Finally, proper quantitative assessment of errors poses a non-trivial task and re- 
mains a challenge for future work, although first steps in this direction have been 
taken in Lanyon et al. (201 1). In contrast to quantum computing, quantum simula- 
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tion is usually not interested in obtaining the many -body wave function in a faithful 
way, but rather aims at an accurate prediction of low order correlation functions, as 
is relevant, for example, for phase diagrams in equilibrium physics. Thus it is gener- 
ally argued that quantum simulation is more robust against errors and imperfections 
than quantum computing, and from an experimental point of view the realization of 
a large-scale quantum simulator is expected to be a more realistic short-term goal 
than building a fault-tolerant quantum computer. However, one of the outstanding 
problems is to investigate the role of errors in an interplay between theory and ex- 
periment. Along a similar line, questions of validation and verification of quantum 
simulators need to be addressed in these future studies. In the context of digital 
quantum simulation, the good news is that - if the gate fidelities and system sizes 
can be further increased - the gate-based approach can incorporate quantum error 
correction protocols. These might prove essential for fault- tolerant quantum simu- 
lation, in particular for future large-scale quantum simulations of complex many- 
body models. 

In the field of quantum information processing, it is one of the grand challenges 
and visions to build in the laboratory a quantum device which performs tasks not 
achievable on a classical level. A next generation quantum simulation experiment 
involving (experimentally proven) large-scale entanglement may be the first labo- 
ratory demonstration that fulfills this promise in a convincing way. This would be 
an exciting and big step forward towards the realization of Feynman's 30-years-old 
dream of building a programmable quantum simulator, which might not only pro- 
vide us with answers to long-standing open questions, but also allow us to explore 
new realms of physics, such as many-body quantum dynamics beyond thermody- 
namic equilibrium. 



5 Acknowledgments 



We acknowledge support by the Austrian Science Fund (FWF) through SFB FO- 
QUS and the START grant Y 58 1-N16 (S. D.), the European Commission (AQUTE, 
NAMEQUAM), the Institut fuer Quanteninformation GmbH and the DARPA OLE 
program. M. M. acknowledges support by the CAM research consortium QUITEMAD 
S2009-ESP-1594, European CommissionPICC:FP7 2007-2013, GrantNo. 249958, 
and the Spanish MICINN grant FIS2009- 10061. 



References 

Abrams, D. S., Lloyd, S., 1997. Simulation of many-body Fermi systems on a uni- 
versal quantum computer. Phys. Rev. Lett. 79 (13), 2586. 
Aguado, M., Brennen, G. K., Verstraete, R, Cirac, J. I., 2008. Creation, manipula- 



67 



tion, and detection of Abelian and non-Abelian anyons in optical lattices. Phys. 
Rev. Lett. 101,260501. 

Alicea, J., Oreg, Y., Refael, G., von Oppen, R, Fisher, M. R A., 201 1. Non-Abelian 
statistics and topological quantum information processing in Id wire networks. 
NaturePhys. 7, 412417. 

Altland, A., Simons, B., 2010. Condensed Matter Field Theory. Cambridge Uni- 
versity Press. 

Altland, A., Zirnbauer, M. R., 1997. Nonstandard symmetry classes in mesoscopic 
normal- superconducting hybrid structures. Phys. Rev. B 55, 1142. 

Altman, E., Auerbach, A., 2002. Plaquette boson-fermion model of cuprates. Phys. 
Rev. B 65, 104508. 

Anderson, P. W., 1987. The resonating valence bond state in La2Cu04 and super- 
conductivity. Science 235, 1196. 

Anderson, S. E., Younge, K. C, Raithel, G., 2011. Trapping Rydberg atoms in an 
optical lattice. Phys. Rev. Lett. 107, 263001. 

Aspect, A., Arimondo, E., Kaiser, R., Vansteenkiste, N., Cohen-Tannoudji, C, 
1988. Laser cooling below the one-photon recoil energy by velocity- selective 
coherent population trapping. Phys. Rev. Lett. 61, 826. 

Bacon, D., Childs, A. M., Chuang, L L., Kempe, J., Leung, D. W., Zhou, X., 2001. 
Universal simulation of Markovian quantum dynamics. Phys. Rev. A 64 (6), 
062302. 

Bakr, W. S., Gillen, J. L, Peng, A., Foiling, S., Greiner, M., 2010a. A quantum 
gas microscope for detecting single atoms in a Hubbard-regime optical lattice. 
Nature 462, 74. 

Bakr, W. S., Peng, A., Tai, M. E., Ma, R., Simon, J., Gillen, J. L, Foiling, S., PoUet, 
L., Greiner, M., 2010b. Probing the superfluid - Mott insulator transition at the 
single-atom level. Science 329 (5991), 547-550. 

Baranov, M. A. Dalmonte, M., Pupillo, G., ZoUer, P., 2012. Condensed matter the- 
ory of dipolar quantum gases. Chemical Reviews, in press. 

Bardyn, C. E., Baranov, M. A., Rico, E., Imamoglu, A., ZoUer, P., Diehl, S., 2012. 
Majorana modes in driven-dissipative atomic superfluids with zero Chem num- 
ber. arXiv: 1201.21 12. 

Barreiro, J., Miiller, M., Schindler, P., Nigg, D., Monz, T., Chwalla, M., Hennrich, 
M., Roos, C. F, ZoUer, P., Blatt, R., 2011. An open-system quantum simulator 
with trapped ions. Nature 470, 486. 

Baumgartner, B., Narnhofer, H., Thirring, W., 2007. Analysis of quantum semi- 
groups with GKS - Lindblad generators: L simple generator. J. Phys. A 41, 
065201. 

Beige, A., Braun, D., Tregenna, B., Knight, P. L., 2000. Quantum computing using 
dissipation to remain in a decoherence-free subspace. Phys. Rev. Lett. 85, 1762- 
1765. 

Benhelm, J., Kirchmair, G., Roos, C. R, Blatt, R., 2008. Towards fault-tolerant 
quantum computing with trapped ions. Nature Phys. 4, 463^66. 

Bermudez, A., Porras, D., Martin-Delgado, M. A., 2009. Competing many-body 
interactions in systems of trapped ions. Phys. Rev. A 79 (6), 060303. 



68 



Berry, D. W., Ahokas, G., Cleve, R., Sanders, B. C, 2007. Efficient quantum al- 
gorithms for simulating sparse Hamiltonians. Communications in Mathematical 
Physics 270, 359. 

Berry, M. V., 1984. Quantal phase factors accompanying adiabatic changes. Proc. 
R. Soc. A 392, 4557. 

Blakestad, R. B., Ospelkaus, C, VanDevender, A. P., Amini, J. M., Britton, J., 
Leibfried, D., Wineland, D. J., 2009. High-fidelity transport of trapped-ion qubits 
through an x-junction trap array. Phys. Rev. Lett. 102 (15), 153002. 

Blatt, R., Wineland, D., 2008. Entangled states of trapped atomic ions. Nature 453, 
1008-1015. 

Bloch, I., Dalibard, J., Nascimbene, S., 2012. Quantum simulations with ultracold 
quantum gases, submitted to Nature Physics. 

Bloch, I., Dalibard, J., Zwerger, W., 2008. Many-body physics with ultracold gases. 
Rev. Mod. Phys. 80, 885. 

Bolognani, S., Ticozzi, R, 2010. Engineering stable discrete-time quantum dynam- 
ics via a canonical QR decomposition. IEEE Trans. Aut. Cont. 55, 2721. 

Bombin, H., Martin-Delgado, M. A., 2006. Topological quantum distillation. Phys. 
Rev. Lett. 97(18), 180501. 

Bombin, H., Martin-Delgado, M. A., 2007. Optimal resources for topological two- 
dimensional stabilizer codes: Comparative study. Phys. Rev. A 76 (1), 012305. 

Brachmann, J. F. S., Bakr, W. S., Gillen, J., Peng, A., Greiner, M., 2011. Induc- 
ing vortices in a Bose-einstein condensate using holographically produced light 
beams. arXiv:1102.3373vl. 

Bravyi, S., DiVincenzo, D. P., Loss, D., Terhal, B. M., 2008. Quantum simulation 
of many -body Hamiltonians using perturbation theory with bounded- strength in- 
teractions. Phys. Rev. Lett. 101 (7), 070503. 

Brennen, G. K., Pupillo, G., Rey, A. M., Clark, C. W., Williams, C. J., 2005. Scal- 
able register initialization for quantum computing in an optical lattice. J. Phys. B 
38, 1687. 

Biichler, H. P., Micheli, A., ZoUer, P., 2007. Three-body interactions with cold polar 
molecules. Nature Phys. 3 (10), 726. 

Buluta, I., Nori, R, 2009. Quantum simulators. Science 326 (5949), 108. 

Calabrese, P., Cardy, J., 2006. Time dependence of correlation functions following 
a quantum quench. Phys. Rev. Lett. 96, 136801. 

Calderbank, A. R., Shor, P. W., 1996. Good quantum error-correcting codes exist. 
Phys. Rev. A 54 (2), 1098-1 105. 

CaroUo, A., Fuentes-Guridi, I., Santos, M. R, Vedral, V., 2003. Geometric phase in 
open systems. Phys. Rev. Lett. 90, 160402. 

Carr, L. D., Ye, J. (Eds.), 2009. Focus on Cold and Ultracold Molecules, New J. 
Phys. Vol. 11. Institute of Physics. 

Ceperley, D. M., 1995. Path integrals in the theory of condensed helium. Rev. Mod. 
Phys. 67, 279-355. 

Chakravarty, S., Ingold, G.-L., Kivelson, S., Luther, A., 1986. Onset of global phase 
coherence in Josephson-junction arrays: A dissipative phase transition. Phys. 
Rev. Lett. 56, 2303-2306. 



69 



Chakravarty, S., Kivelson, S., Zimanyi, G. T., Halperin, B. L, 1987. Effect of quasi- 

particle tunneling on quantum-phase fluctuations and the onset of superconduc- 
tivity in granular films. Phys. Rev. B 35, 7256-7259. 
Chin, J. K., Miller, D. E., Liu, Y., Stan, C, Setiawan, C, Sanner, C, Xu, K., Ket- 

terle, W., 2006. Evidence for superfluidity of ultracold fermions in an optical 

lattice. Nature 443, 961-964. 
Cirac, J. I., ZoUer, P., 2012. Quantum simulation - quo vadis? submitted to Nature 

Physics. 
Clark, R. J., Linb, T., Brown, K. R., Chuang, I. L., 2009. A two-dimensional lattice 

ion trap for quantum simulation. J. Appl. Phys. 105 (1), 013114. 
Clarke, J., Wilhelm, F. K., 2008. Superconducting quantum bits. Nature 453 (7198), 

1031. 
Cramer, M., Dawson, C. M., Eisert, J., Osborne, T. J., 2008. Exact relaxation in a 

class of nonequilibrium quantum lattice systems. Phys. Rev. Lett. 100, 030602. 
Cubel, T., Teo, B. K., Malinovsky, V. S., Guest, J. R., Reinhard, A., Knufli'man, B., 

Berman, P. R., Raithel, G., 2005. Coherent population transfer of ground-state 

atoms into Rydberg states. Phys. Rev. A 72, 023405. 
Daley, A. J., 2011. Quantum computing and quantum simulation with group-ii 

atoms. Quant. Inf. Process. 10, 865-884. 
Daley, A. J., Boyd, M. M., Ye, J., ZoUer, P., 2008. Quantum computing with 

alkaline-earth-metal atoms. Phys. Rev. Lett. 101, 170504. 
Daley, A. J., Taylor, J. M., Diehl, S., Baranov, M., ZoUer, P., 2009. Atomic three- 
body loss as a dynamical three-body interaction. Phys. Rev. Lett. 102, 040402. 
Dalibard, J., Gerbier, E, Juzeliunas, G., Ohberg, P., 2011. Colloquium: Artificial 

gauge potentials for neutral atoms. Rev. Mod. Phys. 83, 1523-1543. 
Dalla Torre, E. G., Demler, E., Giamarchi, T, Altman, E., 2010. Quantum critical 

states and phase transitions in the presence of non-equilibrium noise. Nature 

Physics 6, 806. 
Dalla Torre, E. G., Demler, E., Giamarchi, T, Altman, E., 2011. Dynamics and 

universality in noise driven dissipative systems. arXiv: 11 10.3678. 
Dennis, E., Kitaev, A., Landahl, A., J, P., 2002. Topological quantum memory. J. 

Math. Phys. 43 (9), 4452. 
Diehl, S., Baranov, M., Daley, A. J., ZoUer, P., 2010a. Observability of quantum 

criticality and a continuous supersolid in atomic gases. Phys. Rev. Lett. 104, 

165301. 
Diehl, S., Micheh, A., Kantian, A., Kraus, B., Buchler, H. P, ZoUer, P, 2008. Quan- 
tum states and phases in driven open quantum systems with cold atoms. Nature 

Phys. 4(11), 878. 
Diehl, S., Rico, E., Baranov, M. A., ZoUer, P., 2011. Topology by dissipation in 

atomic quantum wires. Nature Phys. 7, 971977. 
Diehl, S., Tomadin, A., Micheli, A., Fazio, R., ZoUer, P., 2010b. Dynamical phase 

transitions and instabilities in open atomic many-body systems. Phys. Rev. Lett. 

105, 015702. 
Diehl, S., Yi, W., Daley, A. J., ZoUer, P., 2010c. Dissipation-induced J-wave pairing 

of fermionic atoms in an optical lattice. Phys. Rev. Lett. 105, 227001. 



70 



Diir, W., Bremner, M. J., Briegel, H. J., 2008. Quantum simulation of interacting 
high-dimensional systems: The influence of noise. Phys. Rev. A 78 (5), 052325. 

Eisert, J., Prosen, T., 2010. Noise-driven quantum criticality. arXiv: 1012.5013. 

Esslinger, T., 2010. Fermi-Hubbard physics with atoms in an optical lattice. Ann. 
Rev. of Cond. Mat. Phys. 1, 129-152. 

Fisher, M. P A., Weichman, P B., Grinstein, G., Fisher, D. S., 1989. Boson local- 
ization and the superfluid-insulator transition. Phys. Rev. B 40, 546-570. 

Fleischhauer, M., Imamoglu, A., Marangos, J. P., 2005. Electromagnetically in- 
duced transparency: Optics in coherent media. Rev. Mod. Phys. 77, 633. 

Foiling, S., Trotzky, S., Cheinet, P, Feld, M., Saers, R., Widera, A., Miiller, T, 
Bloch, I., 2007. Direct observation of second-order atom tunnelling. Nature 448, 
1029. 

Friedenauer, A., Schmitz, H., Glueckert, J. T., Porras, D., Schaetz, T., 2008. Simu- 
lating a quantum magnet with trapped ions. Nature Phys. 4, 757-761. 

Gaetan, A., Miroshnychenko, Y., Wilk, T, Chotia, A., Viteau, M., Comparat, D., 
Pillet, P., Browaeys, A., Grangier, P., 2009. Observation of collective excitation 
of two individual atoms in the Rydberg blockade regime. Nature Phys. 5, 115. 

Gallagher, T. F, 1994. Rydberg Atoms. Cambridge University Press. 

Gardiner, C. W., ZoUer, P., 1999. Quantum Noise. Springer Verlag, Berlin. 

Gasenzer, T., Berges, J., Schmidt, M. G., Seco, M., 2005. Nonperturbative dynam- 
ical many-body theory of a Bose-Einstein condensate. Phys. Rev. A 72, 063604. 

Georges, A., Kotliar, G., Krauth, W., Rozenberg, M. J., 1996. Dynamical mean- 
field theory of strongly correlated fermion systems and the limit of infinite di- 
mensions. Rev. Mod. Phys. 68, 13-125. 

Gilchrist, A., Langford, N. K., Nielsen, M. A., 2005. Distance measures to compare 
real and ideal quantum processes. Phys. Rev. A 71, 062310. 

Gorshkov, A. V., Hermele, M., Gurarie, V., Xu, C, Julienne, P. S., Ye, J., ZoUer, P., 
Demler, E., Lukin, M. D., Rey, A. M., 2010. Two-orbital SU(N) magnetism with 
ultracold alkaline-earth atoms. Nat. Phys. 6, 289. 

Greiner, M., Mandel, O., Esslinger, T., Hansch, T. W., Bloch, I., 2001. Quantum 
phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms. 
Nature 415 (6867), 39. 

Greiner, M., Mandel, O., Esslinger, T., Hansch, T. W., Bloch, I., 2002a. Quantum 
phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms. 
Nature 415, 39-44. 

Greiner, M., Mandel, O., Hansch, T. W., Bloch, I., 2002b. Collapse and revival of 
the matter wave field of a Bose-Einstein condensate. Nature 419, 51-54. 

Griessner, A., Daley, A. J., Clark, S. R., Jaksch, D., ZoUer, P., 2006. Dark-state 
cooling of atoms by superfluid immersion. Phys. Rev. Lett. 97, 220403. 

Gros, C, 1988. Superconductivity in correlated wave functions. Phys. Rev. B 38, 
931. 

Haffner, H., Roos, C. F, Blatt, R., 2008. Quantum computing with trapped ions. 
Phys. Rep. 469 (4), 155. 

Hallberg, K. A., 2006. New trends in density matrix renormalization. Adv. in Phys. 
55, 477-526. 



71 



Hammond, B. L., Lester, W. A., Reynolds, P. J., 1994. Monte Carlo methods in Ab 
Initio quantum chemistry. World Scientific. 

Hanson, R., Kouwenhoven, L. R, Petta, J. R., Tarucha, S., Vandersypen, L. M. K., 
2007. Spins in few-electron quantum dots. Rev. Mod. Phys. 79 (4), 1217-1265. 

Hasan, M., Kane, C, 2010. Colloquium: Topological insulators. Rev. Mod. Phys. 
82, 3045-3067. 

Hatsugai, Y., 1993. Chern number and edge states in the integer quantum Hall ef- 
fect. Phys. Rev. Lett. 71, 3697-3700. 

Heidemann, R., Raitzsch, U., Bendkowsky, V., Butscher, B., Low, R., Santos, L., 
Pfau, T., 2007. Evidence for coherent collective Rydberg excitation in the strong 
blockade regime. Phys. Rev. Lett. 99, 163601. 

Hensinger, W. K., Olmschenk, S., Stick, D., Hucul, D., Yeo, M., Acton, M., Deslau- 
riers, L., Monroe, C, Rabchuk, J., 2006. T-junction ion trap array for two- 
dimensional ion shuttling, storage, and manipulation. Appl. Phys. Lett. 88 (3), 
034101. 

Hermele, M., Fisher, M. P. A., Balents, L., 2004. Pyrochlore photons: The u(l) spin 
liquid in a 5 = ^ three-dimensional frustrated magnet. Phys. Rev. B 69, 064404. 

Hoening, M., Moos, M., Fleischhauer, M., 2010. Critical exponents of flux- 
equilibrium phase transitions in fermionic lattice models. arXiv: 1 108.2263. 

Hofferberth, S., Lesanovsky, L, Fischer, B., Schumm, T., Schmiedmayer, J., 2007. 
Non-equilibrium coherence dynamics in one-dimensional Bose gases. Nature 
449, 324-327. 

Hofstetter, W., Cirac, J. L, ZoUer, P, Demler, E., Lukin, M. D., 2002. High- 
temperature superfluidity of fermionic atoms in optical lattices. Phys. Rev. Lett. 
89, 220407. 

Home, J. P., Hanneke, D., lost, J. D., Amini, J. M., Leibfried, D., Wineland, D. J., 
2009. Complete methods set for scalable ion trap quantum information process- 
ing. Science 325 (5945), 1227-1230. 

Honer, J., Weimer, H., Pfau, T., Biichler, H. P., 2010. Collective many-body inter- 
action in Rydberg dressed atoms. Phys. Rev. Lett. 105, 160404. 

Isenhower, L., Urban, E., Zhang, X. L., Gill, A. T., Henage, T., Johnson, T. A., 
Walker, T. G., Saffman, M., 2010. Demonstration of a neutral atom controUed- 
NOT quantum gate. Phys. Rev. Lett. 104 (1), 010503. 

Jaksch, D., Bruder, C, Cirac, J. L, Gardiner, C. W., ZoUer, P., 1998a. Cold bosonic 
atoms in optical lattices. Phys. Rev. Lett. 81, 3108. 

Jaksch, D., Bruder, C, Cirac, J. I., Gardiner, C. W., ZoUer, P., 1998b. Cold bosonic 
atoms in optical lattices. Phys. Rev. Lett. 81, 3108-3111. 

Jaksch, D., Cirac, J. I., ZoUer, P, Rolston, S. L., Cote, R., Lukin, M. D., 2000. Fast 
quantum gates for neutral atoms. Phys. Rev. Lett. 85 (10), 2208. 

Jaksch, D., ZoUer, P, 2005. The cold atom Hubbard toolbox. Ann. Phys. 315 (1), 
52. 

Jane, E., Vidal, G., Diir, W., ZoUer, P., Cirac, J., 2003. Simulation of quantum dy- 
namics with quantum optical systems. Quantum Information & Computation 3, 
15-37. 

Johanning, M., Vam, A. F, Wunderlich, C, 2009. Quantum simulations with cold 



72 



trapped ions. J. Phys. B 42 (15), 154009. 

Jordens, R., Strohmaier, N., Giinter, K., Moritz, H., Esslinger, T., 2008. Fermionic 
atoms in a three dimensional optical lattice: Observing Fermi surfaces, dynamics, 
and interactions. Nature 455, 204-207. 

Jordens, R., Tarruell, L., Greif, D., Uehlinger, T., Strohmaier, N., Moritz, H., 
Esslinger, T., De Leo, L., KoUath, C, Georges, A., Scarola, V., PoUet, L., 
Burovski, E., Kozik, E., Troyer, M., 2010. Quantitative determination of tem- 
perature in the approach to magnetic order of ultracold fermions in an optical 
lattice. Phys. Rev. Lett. 104, 180401. 

Jozsa, R., 1994. Fidelities for mixed quantum states. J. Mod. Opt. 41, 2315. 

Kamenev, A., Levchenko, A., 2009. Keldysh technique and non-linear sigma- 
model: basic principles and applications. Adv. in Phys. 58, 197. 

Kampf, A., Schon, G., 1987. Quantum effects and the dissipation by quasiparticle 
tunneling in arrays of Josephson junctions. Phys. Rev. B 36, 3651-3660. 

Kantian, A., Dalmonte, M., Diehl, S., Hofstetter, W., ZoUer, P., Daley, A. J., 2009. 
Atomic color superfluid via three-body loss. Phys. Rev. Lett. 103, 240401. 

Kasevich, M., Chu, S., 1992. Laser cooling below a photon recoil with three-level 
atoms. Phys. Rev. Lett. 69, 1741. 

Kasprzak, J., Richard, M., Kundermann, S., Baas, A., Jeambrun, P., Keeling, J. 
M. J., Marchetti, F M., Szymaska, M. H., Andre, R., Staehli, J. L., Savona, V., 
Littlewood, P. B., Deveaud, B., Dang, L. S., 2006. Bose-Einstein condensation 
of exciton polaritons. Nature 443, 409. 

Kassal, I., Whitfield, J. D., Perdomo-Ortiz, A., Yung, M.-H., Aspuru-Guzik, A., 
2011. Simulating chemistry using quantum computers. Annual Review of Phys- 
ical Chemistry 62, 185. 

Kastoryano, M. J., Reiter, R, S0rensen, A. S., 2011. Dissipative preparation of en- 
tanglement in optical cavities. Phys. Rev. Lett. 106, 090502. 

Kinoshita, T., Wenger, T., Weiss, D. S., 2006. A quantum Newton's cradle. Nature 
452, 900-903. 

Kitaev, A. Y., 1997. Quantum error correction with imperfect gates. Russ. Math. 
Surv. 52, 1191. 

Kitaev, A. Y, 2001. Unpaired Majorana fermions in quantum wires. Physics- 
Uspekhi44, 131. 

Kitaev, A. Y, 2003. Fault-tolerant quantum computation by anyons. Ann. Phys. 
303(1), 2. 

Kitaev, A. Y, 2006. Anyons in an exactly solved model and beyond. Annals of 
Physics 321, 2. 

Kitagawa, T., Berg, E., Rudner, M., Demler, E., 2010. Topological characterization 
of periodically driven quantum systems. Phys. Rev. B 82, 235114. 

Kliesch, M., Barthel, T., Gogolin, C, Kastoryano, M., Eisert, J., 2011. Dissipative 
quantum Church-Turing theorem. Phys. Rev. Lett. 107, 120501. 

Kogut, J. B., 1979. An introduction to lattice gauge theory and spin systems. Rev. 
Mod. Phys. 51,659-713. 

Kohl, M., Moritz, H., Stoferle, T, Giinter, K., Esslinger, T., 2005. Fermionic atoms 
in a three dimensional optical lattice: Observing Fermi surfaces, dynamics, and 



73 



interactions. Phys. Rev. Lett. 94, 080403. 
Kollath, C, Lauchli, A. M., Altman, E., 2007. Quench dynamics and nonequilib- 

rium phase diagram of the Bose-Hubbard model. Phys. Rev. Lett. 98, 18060L 
Kraus, B., Buchler, H. P, Diehl, S., Kantian, A., Micheli, A., Zoller, P, 2008. Prepa- 
ration of entangled states by quantum Markov processes. Phys. Rev. A 78 (4), 

042307. 
Kraus, C. V., Diehl, S., Zoller, P., Baranov, M. A., 2012. Probing atomic Majorana 

fermions in optical lattices. arXiv: 1201.3253. 
Krauter, H., Muschik, C. A., Jensen, K., Wasilewski, W., Petersen, J. M., Cirac, 

J. L, Polzik, E. S., 2011. Entanglement generated by dissipation and steady state 

entanglement of two macroscopic objects. Phys. Rev. Lett. 107, 080503. 
Ladd, T. D., Jelezko, R, Laflamme, R., Nakamura, Y., Monroe, C, O'Brien, J. L., 

2010. Quantum computers. Nature 464 (7285), 45-53. 
Lanyon, B. P., Hempel, C, Nigg, D., Miiller, M., Gerritsma, R., Zahringer, R, 

Schindler, P., Barreiro, J. T., Rambach, M., Kirchmair, G., Hennrich, M., Zoller, 

P., Blatt, R., Roos, C. E, 201 1. Universal digital quantum simulation with trapped 

ions. Science 334 (6052), 57. 
Lee, P. J., Anderlini, M., Brown, B. L., Sebby-Strabley, J., Phillips, W. D., Porto, 

J. v., 2007. Sublattice addressing and spin-dependent motion of atoms in a 

double-well lattice. Phys. Rev. Lett. 99, 020402. 
Leggett, A. J., 2006. Quantum Liquids: Bose Condensation and Cooper Pairing in 

Condensed-Matter Systems. Oxford University Press. 
Lesanovsky, L, 2011. Many-body spin interactions and the ground state of a dense 

Rydberg lattice gas. Phys. Rev. Lett. 106, 025301. 
Levin, M., Wen, X.-G., 2005a. Colloquium: Photons and electrons as emergent 

phenomena. Rev. Mod. Phys. 77, 871-879. 
Levin, M. A., Wen, X.-G., 2005b. String-net condensation: A physical mechanism 

for topological phases. Phys. Rev. B 71, 0451 10. 
Lewenstein, M., Sanpera, A., Ahufinger, V., Damski, B., De, A. S., Sen, U., 2006. 

Ultracold atomic gases in optical lattices: mimicking condensed matter physics 

and beyond. Advances in Physics 56, 135. 
Lidar, D. A., Chuang, L L., Whaley, K. B., 1998. Decoherence-free subspaces for 

quantum computation. Phys. Rev. Lett. 81, 2594-2597. 
Lindner, N. H., Refael, G., Galitski, V., 201 1 . Eloquet topological insulator in semi- 
conductor quantum wells. Nature Phys. 7, 490495. 
Lloyd, S., 1995. Almost any quantum logic gate is universal. Phys. Rev. Lett. 75, 

346-349. 
Lloyd, S., 1996. Universal quantum simulators. Science 273 (5278), 1073. 
Lloyd, S., Viola, L., 2001. Engineering quantum dynamics. Phys. Rev. A 65 (1), 

010101. 
Lukin, M. D., Eleischhauer, M., Cote, R., Duan, L. M., Jaksch, D., Cirac, J. L, 

Zoller, P., 2001. Dipole blockade and quantum information processing in meso- 

scopic atomic ensembles. Phys. Rev. Lett. 87 (3), 037901. 
Makela, H., Watanabe, G., 2011. Dissipation-induced squeezing. 

arXiv:1101.4845vl. 



74 



Marcos, D., Tomadin, A., Diehl, S., Rabl, P., 2012. Photon condensation in circuit 
QED by engineered dissipation. arxiv:1201.3388vl. 

Mariantoni, M., Wang, H., Bialczak, R. C, Lenander, M., Lucero, E., Neeley, M., 
OConnell, A. D., Sank, D., Weides, M., Wenner, J., Yamamoto, T., Yin, Y, Zhao, 
J., Martinis, J. M., Cleland, A. N., 2011. Photon shell game in three-resonator 
circuit quantum electrodynamics. Nature Physics 7, 287293. 

Mitra, A., Takei, S., Kim, Y B., Millis, A. J., 2006. Nonequilibrium quantum criti- 
cality in open electronic systems. Phys. Rev. Lett. 97, 236808. 

Moessner, R., Sondhi, S. L., 2001. Resonating valence bond phase in the triangular 
lattice quantum dimer model. Phys. Rev. Lett. 86 (9), 1881-1884. 

Mohapatra, A. K., Jackson, T. R., Adams, C. S., 2007. Coherent optical detection 
of highly excited Rydberg states using electromagnetically induced transparency. 
Phys. Rev. Lett. 98, 113003. 

M0ller, D., Madsen, L. B., M0lmer, K., 2008. Quantum gates and multiparticle 
entanglement by Rydberg excitation blockade and adiabatic passage. Phys. Rev. 
Lett. 100, 170504. 

M0lmer, K., S0rensen, A., 1999. Multiparticle entanglement of hot trapped ions. 
Phys. Rev. Lett. 82 (9), 1835. 

Monz, T., Schindler, P, Barreiro, J. T, Chwalla, M., Nigg, D., Coish, W. A., Har- 
lander, M., Hansel, W., Hennrich, M., Blatt, R., 2011. 14-qubit entanglement: 
Creation and decoherence. Phys. Rev. Lett. 106 (13), 130605. 

Moskalenko, S. A., Snoke, D., 2000. Bose-Einstein Condensation of Excitons and 
Biexcitons. Cambridge Univ. Press. 

Motrunich, O. I., Senthil, T., 2002. Exotic order in simple models of bosonic sys- 
tems. Phys. Rev. Lett. 89, 277004. 

Miiller, M., Hammerer, K., Zhou, Y, Roos, C. R, ZoUer, P., 2011. Simulating open 
quantum systems: from many-body interactions to stabilizer pumping. New J. 
Phys. 13, 085007. 

MuUer, M., Lesanovsky, I., Weimer, H., Buchler, H. P., ZoUer, P., 2009. Meso- 
scopic Rydberg gate based on electromagnetically induced transparency. Phys. 
Rev. Lett. 102, 170502. 

Muschik, C. A., Polzik, E. S., Cirac, J. I., 2011. Dissipatively driven entanglement 
of two macroscopic atomic ensembles. Phys. Rev. A 83, 052312. 

Nayak, C, Simon, S. H., Stem, A., Freedman, M., Das Sarma, S., 2008. Non- 
Abelian anyons and topological quantum computation. Rev. Mod. Phys. 80 (3), 
1083. 

Nelson, K. D., Li, X., Weiss, D. S., 2007. Imaging single atoms in a three- 
dimensional array. Nature Phys. 3, 556. 

Nielsen, M. A., Chuang, I. L., 2000. Quantum Computation and Quantum Infor- 
mation. Cambridge University Press. 

O'Brien, J. L., 2007. Optical quantum computing. Science 318, 1567-1570. 

Olmos, B., Gonzalez-Ferez, R., Lesanovsky, I., 2009. Fermionic collective excita- 
tions in a lattice gas of Rydberg atoms. Phys. Rev. Lett. 103, 185302. 

Ortiz, G., Gubernatis, J. E., Knill, E., Laflamme, R., 2001. Quantum algorithms for 
fermionic simulations. Phys. Rev. A 64 (2), 022319. 



75 



Pachos, J., Zanardi, P., Rasetti, M., 1999. Non-Abelian Berry connections for quan- 
tum computation. Phys. Rev. A 61, 010305. 

Paramekanti, A., Randeria, M., Trivedi, N., 2004. High-Tc superconductors? A 
variational theory of the superconducting state. Phys. Rev. B 70, 054504. 

Parr, R. G., Yang, W., 1989. Density-Functional Theory of Atoms and Molecules. 
Oxford University Press. 

Pastawski, R, Clemente, L., Cirac, J. I., 2011. Quantum memories based on engi- 
neered dissipation. Phys. Rev. A 83, 012304. 

Pichler, H., Daley, A. J., ZoUer, P., 2010. Nonequilibrium dynamics of bosonic 
atoms in optical lattices: Decoherence of many-body states due to spontaneous 
emission. Phys. Rev. A 82, 063605. 

Pohl, T., Demler, E., Lukin, M. D., 2010. Dynamical crystallization in the dipole 
blockade of ultracold atoms. Phys. Rev. Lett. 104, 043002. 

Porras, D., Cirac, J. I., 2004. Effective quantum spin systems with trapped ions. 
Phys. Rev. Lett. 92 (20), 207901. 

Porto, J. v., 2008. Improving correlations despite particle loss. Science 320, 1300. 

Prokof'ev, N., Svistunov, B., Tupitsyn, I., 1998. Exact, complete, and universal 
continuous-time worldline Monte Carlo approach to the statistics of discrete 
quantum systems. JETP 87, 310. 

Prosen, T., 2008. Third quantization: A general method to solve master equations 
for quadratic open Fermi systems. New J. Phys. 10, 043026. 

Prosen, T., 2010. Spectral theorem for the Lindblad equation for quadratic open 
fermionic systems. J. Stat. Mech., P07020. 

Pupillo, G., Rey, A. M., Brennen, G. K., Williams, C. J., Clark, C. W., 2004. Scal- 
able quantum computation in systems with Bose-Hubbard dynamics. J. Mod. 
Opt. 51,2395. 

Qi, X.-L., Zhang, S.-C, 2011. Topological insulators and superconductors. Rev. 
Mod. Phys. 83, 1057-1110. 

Raitzsch, U., Bendkowsky, V., Heidemann, R., Butscher, B., Low, R., Pfau, T, 
2008. Echo experiments in a strongly interacting Rydberg gas. Phys. Rev. Lett. 
100, 013002. 

Read, N., Green, D., 2000. Paired states of fermions in two dimensions with break- 
ing of parity and time-reversal symmetries and the fractional quantum Hall effect. 
Phys. Rev. B 61, 10267-10297. 

Reetz-Lamour, M., Amthor, T., Deiglmayr, J., Weidemiiller, M., 2008. Rabi oscil- 
lations and excitation trapping in the coherent excitation of a mesoscopic frozen 
Rydberg gas. Phys. Rev. Lett. 100, 253001. 

Reichenbach, I., Deutsch, I. H., 2007. Sideband cooling while preserving coher- 
ences in the nuclear spin state in group-ii-like atoms. Phys. Rev. Lett. 99, 123001 . 

Riebe, M., Monz, T., Kim, K., Villar, A. S., Schindler, P., Chwalla, M., Hennrich, 
M., Blatt, R., 2008. Deterministic entanglement swapping with an ion-trap quan- 
tum computer. Nat. Phys. 4, 839. 

Rigol, M., Dunjko, V., Olshanii, M., 2008. Thermalization and its mechanism for 
generic isolated quantum systems. Nature 452, 854-858. 

Rokhsar, D. S., Kivelson, S. A., 1988. Superconductivity and the quantum hard- 



76 



core dimer gas. Phys. Rev. Lett. 61, 2376-2379. 

Roncaglia, M., Rizzi, M., Cirac, J. I., 2010. PfafRan state generation by strong 
three-body dissipation. Phys. Rev. Lett. 104, 096803. 

Roos, C. R, 2008. Ion trap quantum gates with amplitude-modulated laser beams. 
New J. Phys. 10, 013002. 

Rudner, M. S., Levitov, L. S., 2009. Topological transition in a non-hermitian quan- 
tum walk. Phys. Rev. Lett. 102, 065703. 

Ryu, S., Schnyder, A., Furusaki, A., Ludwig, A. W. W., 2010. Topological insula- 
tors and superconductors: Ten-fold way and dimensional hierarchy. New J. Phys. 
12, 065010. 

Sachdev, S., 1999. Quantum Phase Transitions. Cambridge University Press. 

Sadler, L. E., Higbie, J. M., Leshe, S. R., Vengalattore, M., Stamper-Kurn, D. M., 
2006. Spontaneous symmetry breaking in a quenched ferromagnetic spinor 
Bose-Einstein condensate. Nature 443, 312-315. 

Saffman, M., Walker, T. G., M0lmer, K., 2010. Quantum information with Rydberg 
atoms. Rev. Mod. Phys. 82 (3), 2313. 

Schachenmayer, J., Lesanovsky, I., Micheli, A., Daley, A. J., 2010. Dynamical crys- 
tal creation with polar molecules or Rydberg atoms in optical lattices. New J. 
Phys. 12, 103044. 

Schindler, P., Barreiro, J. T, Monz, T., Nebendahl, V., Nigg, D., Chwalla, M., Hen- 
nrich, M., Blatt, R., 201 1 . Experimental repetitive quantum error correction. Sci- 
ence 332, 1059-1061. 

Schmid, A., 1983. Diffusion and localization in a dissipative quantum system. Phys. 
Rev. Lett. 51, 1506-1509. 

Schmied, R., Wesenberg, J. H., Leibfried, D., 2009. Optimal surface-electrode trap 
lattices for quantum simulation with trapped ions. Phys. Rev. Lett. 102, 233002. 

Schneider, C, Porras, D., Schatz, T., 2012. Experimental quantum simulations of 
many-body physics with trapped ions. Reports on Progress in Physics 75 (2), 
024401. 

Schneider, U., Hackermiiller, L., Will, S., Best, T., Bloch, I., Costi, T. A., Helmes, 
R. W., Rasch, D., Rosch, A., 2008. Metallic and insulating phases of repulsively 
interacting fermions in a 3d optical lattice. Science 322 (5907), 1520-1525. 

SchoUwock, U., 2005. The density-matrix renormalization group. Rev. Mod. Phys. 
77,259-315. 

Shavitt, L, Bartlett, R. J., 2009. Many-Body Methods in Chemistry and Physics: 
MBPT and Coupled-Cluster Theory. Cambridge University Press. 

Sherson, J. R, Weitenberg, C, Endres, Manuel Cheneau, M., Bloch, I., Kuhr, S., 
2010. Single-atom-resolved fluorescence imaging of an atomic Mott insulator. 
Nature 467, 68. 

Simon, B., 1983. Holonomy, the quantum adiabatic theorem, and Berry's phase. 
Phys. Rev. Lett. 51, 2167-2170. 

Singer, K., Reetz-Lamour, M., Amthor, T., Marcassa, L. G., Weidemiiller, M., 2004. 
Suppression of excitation and spectral broadening induced by interactions in a 
cold gas of Rydberg atoms. Phys. Rev. Lett. 93, 163001. 

S0rensen, A., M0lmer, K., 1999. Spin-spin interaction and spin squeezing in an 



77 



optical lattice. Phys. Rev. Lett. 83, 2274-2277. 

Steane, A. M., 1996. Error correcting codes in quantum theory. Phys. Rev. Lett. 
77 (5), 793. 

Sun, B., Robicheaux, E, 2008. Numerical study of two-body correlation in a Id 
lattice with perfect blockade. New J. Phys. 10, 045032. 

Suzuki, M., 1992. General theory of higher-order decomposition of exponential 
operators and symplectic integrators. Phys. Lett. A 165, 387. 

Syassen, N., Bauer, D. M., Lettner, M., Volz, T, Dietze, D., Garcia-RipoU, J. J., 
Cirac, J. L, Rempe, G., Diirr, S., 2008. Strong dissipation inhibits losses and 
induces correlations in cold molecular gases. Science 320, 1329-1331. 

Tomadin, A., Diehl, S., ZoUer, P., 201 1. Nonequilibrium phase diagram of a driven 
and dissipative many-body system. Phys. Rev. A 83, 01361 1. 

Tong, D., Earooqi, S. M., Stanojevic, J., Krishnan, S., Zhang, Y. P., Cote, R., Eyler, 
E. E., Gould, P. L., 2004. Local blockade of Rydberg excitation in an ultracold 
gas. Phys. Rev. Lett. 93, 063001. 

Trotzky, S., Chen, Y.-A., Elesch, A., McCuUoch, L P., SchoUwock, U., Eisert, 
J., Bloch, L, 2011. Probing the relaxation towards equilibrium in an isolated 
strongly correlated Id Bose gas. arXiv: 1101.2659. 

Urban, E., Johnson, T. A., Henage, T., Isenhower, L., Yavuz, D. D., Walker, T. G., 
Saffman, M., 2009. Observation of Rydberg blockade between two atoms. Na- 
ture Phys. 5, 110. 

Vandersypen, L. M. K., Chuang, L L., 2004. NMR techniques for quantum control 
and computation. Rev. Mod. Phys. 76 (4), 1037. 

Verstraete, E, Wolf, M. M., Cirac, J. L, 2009. Quantum computation and quantum- 
state engineering driven by dissipation. Nature Phys. 5 (9), 633. 

Vidal, J., Dusuel, S., Schmidt, K. P., 2009. Low-energy effective theory of the toric 
code model in a parallel magnetic field. Phys. Rev. B 79, 033109. 

Viteau, M., Bason, M. G., Radogostowicz, J., Malossi, N., Ciampini, D., Morsch, 
O., Arimondo, E., 2011. Rydberg excitations in Bose-Einstein condensates 
in quasi-one-dimensional potentials and optical lattices. Phys. Rev. Lett. 107, 
060402. 

Vogt, T, Viteau, M., Zhao, J., Chotia, A., Comparat, D., Pillet, P., 2006. Dipole 
blockade at Eorster resonances in high resolution laser excitation of Rydberg 
states of cesium atoms. Phys. Rev. Lett. 97, 083003. 

Walker, T. G., Saffman, M., 2012. Entanglement of two atoms using Rydberg block- 
ade. Advances of Atomic, Molecular and Optical Physics. 

Weimer, H., Low, R., Pfau, T, Biichler, H. P., 2008. Quantum critical behavior in 
strongly interacting Rydberg gases. Phys. Rev. Lett. 101, 250601. 

Weimer, H., Miiller, M., Biichler, H. P., Lesanovsky, L, 201 1. Digital quantum sim- 
ulation with Rydberg atoms. Quantum Inf. Process. 10, 885. 

Weimer, H., Miiller, M., Lesanovsky, L, ZoUer, P., Biichler, H. P., 2010. A Rydberg 
quantum simulator. Nature Phys. 6, 382-388. 

Weitenberg, C, Endres, M., Sherson, J. E, Cheneau, M., Schausz, P., Eukuhara, 
T., Bloch, L, Kuhr, S., 2011. Single-spin addressing in an atomic Mott insulator. 
Nature471 (7338), 319. 



78 



Wen, X.-G., 2004. Quantum Field Theory of Many-Body Systems: From the Origin 
of Sound to an Origin of Light and Electrons. Oxford University Press. 

Whitlock, S., Gerritsma, R., Fernholz, T, Spreeuw, R. J. C, 2009. Two-dimensional 
array of microtraps with atomic shift register on a chip. New J. Phys. 11 (2), 
023021. 

Wilczek, F., Zee, A., 1984. Appearance of gauge structure in simple dynamical 
systems. Phys. Rev. Lett. 52, 2111-2114. 

Wilk, T., Gaetan, A., Evellin, C, Wolters, J., Miroshnychenko, Y., Grangier, P., 
Browaeys, A., 2010. Entanglement of two individual neutral atoms using Ryd- 
berg blockade. Phys. Rev. Lett. 104 (1), 010502. 

Wiseman, H. M., Milbum, G. J., 2009. Quantum Measurement and Control. Cam- 
bridge University Press. 

Wrachtrup, J., Jelezko, R, 2006. Processing quantum information in diamond. J. 
Phys.: Condens. Matter 18, 807. 

Wu, R., Pechen, A., Brif, C, Rabitz, H., 2007. Controllability of open quantum 
systems with Kraus-map dynamics. J. Phys. A: Mathematical and Theoretical 
40 (21), 5681. 

Ye, J., Kimble, H. J., Katori, H., 2008. Quantum state engineering and precision 
metrology using state-insensitive light traps. Science 320, 1734-1738. 

Yi, W., Diehl, S., Daley, A. J., ZoUer, P., 201 1. Driven-dissipative many-body pair- 
ing states for cold fermionic atoms in an optical lattice. arXiv: 1 1 1 1.7053. 

Zhang, F C, Gros, C, Rice, T M., Shiba, H., 1988. A renormalised Hamiltonian 
approach to a resonant valence bond wavefunction. Supercond. Sci. Technol. 1, 
36. 



79 



